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MULTIPLE OBSERVATIONAL CRITERIA: THE CASES OF Z = 6.7 GRB 080913, Z = 8.3 GRB 090423, AND 

SOME SHORT/HARD GRBS 

Bing Zhang 1 , Bin-Bin Zhang 1 , Francisco J. Virgili 1 , En- Wei Liang 2 , D. Alexander Kann 3 , Xue-Feng Wu 4 ' 5 , 
Daniel Proga 1 , Hou-Jun Lv 2 , Kenji Toma 4 , Peter Meszaros 4 ' 6 , David N. Burrows 4 , Peter W. A. Roming 4 , Neil 

Gehrels 7 

Draft version August 6, 2009 

ABSTRACT 

The two high-redshift gamma-ray bursts, GRB 080913 at z = 6.7 and GRB 090423 at z = 8.3, 
recently detected by Swift appear as intrinsically short, hard GRBs. They could have been recognized 
by BATSE as short /hard GRBs should they have occurred at z < 1. In order to address their physi- 
cal origin, we perform a more thorough investigation on two physically distinct types (Type I/II) of 
cosmological GRBs and their observational characteristics. We reiterate the definitions of Type I/II 
GRBs and then review the following observational criteria and their physical motivations: supernova 
association, specific star forming rate of the host galaxy, location offset, duration, hardness, spectral 
lag, statistical correlations, energetics and collimation, afterglow properties, redshift distribution, lu- 
minosity function, and gravitational wave signature. Contrary to the traditional approach of assigning 
the physical category based on the gamma-ray properties (duration, hardness, and spectral lag), we 
take an alternative approach to define the Type I and Type II Gold Samples using several criteria 
that are more directly related to the GRB progenitors (supernova association, host galaxy type, and 
specific star forming rate). We then study the properties of the two Gold Samples and compare them 
with the traditional long/soft and short/hard samples. We find that the Type II Gold Sample reason- 
ably tracks the long/soft population, although it includes several intrinsically short (shorter than Is 
in the rest frame) GRBs. The Type I Gold Sample only has 5 GRBs, 4 of which are not strictly short 
but have extended emission. Other short /hard GRBs detected in the Swift era represent the BATSE 
short/hard sample well, but it is unclear whether all of them belong to Type I. We suggest that some 
(probably even most) high-luminosity short/hard GRBs instead belong to Type II. Based on multiple 
observational criteria, we suggest that GRB 080913 and GRB 090423 are more likely Type II events. 
In general, we acknowledge that it is not always straightforward to discern the physical categories of 
GRBs, and re-emphasize the importance of invoking multiple observational criteria. We cautiously 
propose an operational procedure to infer the physical origin of a given GRB with available multiple 
observational criteria, with various caveats laid out. 

Subject headings: gamma-rays: bursts — gamma rays: observations — gamma rays: theory 



1. INTRODUCTION 

Phenomenologically, gamma-ray bursts (GRBs) have 
been generally classified into the long-duration, soft- 
spectrum class and the short-duration, hard-spectrum 
class in the CGRO/BATSE era based on the bimodal 
distribution of GRBs in the duration-hardness diagram 
(|Kouveliotou et al.lll993l ) 8 . There is no clear boundary 
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8 Several analyses hav e suggested the existence of an inter- 
mediate duration g roup (Mukhcricc ct al. 1998; Horvath 1998; 
Hakkila et al 2000). However, as discussed in the bulk of the text 
below, there is so far no strong indication of the existence of a 
third, physically distinct category of cosmological GRBs based on 



line in this diagram to separate the two populations. 
Traditionally, an observer-frame BATSE-band duration 
Tgo ~ 2 s has been taken as the separation line: bursts 
with Tgo > 2s are "long" and bursts with Tgo < 2s are 
"short". 

The journey was long to uncover the physical ori- 
gins of these two phenomenologically different classes 
of GRBs. The discoveries and the routine observa- 
tions of the broad band afterglows of long GRBs re- 
veal that their host galaxies are typically irregular (in 
a few cases spiral) g alaxies with intense star formation 
(jFruchter et alJl2006h . In a handful of cases these GRBs 
are firmly associated with Type Ib/c supernovae (e.g. 
i Hiorth et all 120031: IStanek et aU 120031 : ICampana et all 
l2006t|Pian et al.|[2006[ ). This strongly suggests that they 
are likely related to deaths of massive stars. Theoreti- 
cally, the "collapsar" model of GRBs has been discussed 
over the years as the standard scenario for long GRBs 
(IWoos lcv 1993 ; |Paczv.ns ki 1998; iMacFadven fc Wooslevl 
ll999HWooslev fc Bloomll2006l ). 

The breakthrough to understand the nature of some 

multiple observational data. So we will focus on the two main 
phenomcnological categories of GRBs in the rest of the paper. 
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short GRBs was made in 2005 after the launch of 
the Swift satellite (jGehrels et all 120041 ). Prompt lo- 
calizations and deep afterglow searches for a hand- 
ful o f short GRBs (iGehrels et alj 120051: IBloom et al 



2006: Fox ct al. 2005; Villascnor ct al. 2005; Hiorth ct al 



2005al iBarthelmv et all l2005al iBerger et all 120051 ) sug- 



gest that some of them are associated with nearby early- 
type galaxies with little star formation. Deep searches 
of associated supernovae from these events all led to 
non-detections (e.g. Kann et al. 2008 and references 
therein, see also Appendix for more references). These 
are in stark contrast to the bursts detected in the pre- 
Swift era (mostly long-duration). On the other hand, the 
observations are consistent with (although not a direct 
proof of) the long-sought progenitor models that invoke 
mergers of two compact stellar objects , leading candi- 
dates being NS-NS and NS-BH systems (lPaczvnskilll986l ; 
lEichler et al.lll989t[Paczvnskill991l;lNaravan et al.lll992f ). 
Although the sample with secure host galaxies is small, 
a general trend in the community is to accept that the 
B ATSE short /hard population bursts are of this compact 
star merger origin 9 . 

The clean dichotomy of the two populations (both 
phenomenological and physical) was soon muddled by 
the detection of a nearby long-durat i on GRB with- 
out S N association (|Gehrels et a l. 2006; Gal- Ya m et al.1 



2006; iFvnbo et all [20061 : 
GRB 060614 has 



IDella Valle elaTl l2006al) . 

100s in the Swift BAT 



(|Barthehny et al. 1 l2005bf ) band, which phenomenologi- 
cally definitely belongs to the long duration category. 
On the other hand, the light curve is characterized by 
a short/hard spike (with a duration ~ 5 s) followed by 
a series of soft gamma-ray pulses. The spectral lag at 
the short/hard spike is negligi bly small, a comm on fea- 
ture of the short /hard GRBs (jGehrels et al.ll2006f) . Very 
stringent upper limits on the radia tion flux from an un- 
derlying SN have been established dGal-Yam et al.l 120061 : 



_JG_ 

IFvnbo et al.ll2006t IDella Valle et al.ll2006af ). These facts 



are consistent with the compact star merger scenario. 
More interestingly, this burst looks like a more ener- 
getic version of GRB 050724, the "sm oking-gun" burst 
of the compact star merge r population ( Barthelmv et al.l 



l2005aHBerger et al.ll^rjll ). IZhang et al l (|2007bl ) showed 



that if one applies the E p oc E^f so relation (jAmati et al.1 
I2002t iLiang et al.l l2004f) to GRB 060614 and makes it 
as energetic as GRB 050724, the pseudo-burst would 
be detected as a marginal short/hard burst by BATSE, 
and would be very similar to GRB 050724 if detected 
by Swift/BAT. In particular, the soft gamma-ray tail 
would appear as the "extended emission" detected in 
some "short/hard" GRBs including GRB 050724. A 
second, much shorter (with Tgo ~ 4 s) burst with- 
out SN associa tion, GRB 060505, was detected around 
the same time ([Fvnbo et al.l 120061 ). However, the phys- 
ical nature of this burst is subject to intense de- 



bate (lOfek et al.||2007: 


Thone et al.ll2008; McBreen et all 


20081: iKann et al.l 2008 


). 



9 It is widely accepted that at least a fraction of short/hard 
GRBs ar e the giant flares of soft gamma-ray r epeaters in nearby 
galaxies (Palmer ct al. 2005; Tanvir ct al. 2005). The observations 
sugges t that the contribu tion fro m such a populati o n is n ot signif- 
icant l INakar et al.l 120061) . but see IChapman et al.1 112003) . We do 
not discuss these bursts in this paper. 



In any case, duration and hardness are not neces- 
sarily reliable indicators of the physical nature of a 
GRB any more. In order to determine whether or 
not a GRB can be associated with a particular phys- 
ical model, one i s forced to appeal to multiple obser- 
vational criteria (|Donaghv et alj I2006D. Prompted by 
the detectio n of GRB 060614, we ([Zhang eFall l2007bl ; 
Zhang 2006) suggested naming the bursts that are con- 
sistent with the massive-star origin and the compact- 
star-merger origin models as Type II and Type I, respec- 
tively 10 , and attempted to invoke a set of multiple obser- 
vational criteria to judge the physical category of a GRB. 
A more dev eloped physic a l cate gorization scheme was 
proposed by IBloom et alj (|2008l ). who also introduced 
SGR giant-flare-like (non-destructive and l ikely repeat- 
ing) e vents. Within the destructive events, IBloom et al.l 
( 2008) agreed that there are two major model types (de- 
generate and non-de generate), wh i ch corr espond to Type 
I and Type II in the lZhang et alj (|2007bf )'s classification 
scheme. Throughout this paper we will adopt the nomen- 
clature of Type I/II to denote the two physically distinct 
categories of cosmological GRB models. 

The recent ly detected two high- z GRBs, GRB 080913 
at z = 6 7 (iGreiner et al.l I2q09al) and GRB 090423 at 
z = 8.3 (jTanvir et all 120091: ISalvaterra et all 120091) in- 
troduce a further complication to the scheme associ- 
ating GRBs with particular theoretical models. Be- 
ing the two GRBs with the highest redshifts as of the 
time of writing, these two bursts each have a redshift- 
corrected duration [Tgo/(l + z)\ shorter than 2 seconds, 
with a hard spectrum typical for short/hard GRBs. 
This naturally raises the interesting qu estion regard- 
ing the progenitor system of the burst (Grciiicr et al.l 
2009aHPerez- Ramirez et al.ll2008l: iBelczvnski et al.ll2008l: 



Tanvir et all 120091 : ISalvaterra et al.ll2009f) . More gener 



ally, it again raises the difficult question regarding how to 
use the observed properties to judge the physical origin 
of a GRB. In this paper, we make some attempts to ad- 
dress this difficult problem. The structure of the paper is 
the following. In §2, we present the observational prop- 
erties of GRB 080913 and GRB 090423, and show that 
if the identical bursts had occurred at z < 1, they could 
have been recognized as short hard GRBs based on their 
observed properties. In §3, we comment on the strengths 
and weaknesses of classifying GRBs based on physically 
motivated criteria. We then reiterate the definitions of 
Type I/II GRBs in §4, and critically review a list of ob- 
servational criteria as well as their physical motivations 
as discussed in the literature. In order to address the pro- 
found questions of whether "Type I = short /hard/short 
lag" and "Type II = long/soft/long lag", in §5 we take 
an alternative approach (from the traditional one) to as- 
sess the problem. Instead of associating a burst with 
a particular physical model (massive star core collapses 
vs. compact star mergers) a priori based on its gamma- 
ray properties (duration, hardness, spectral lag), we use 
several observational properties that are more directly 
relevant to the GRB progenitors to define the Gold Sam- 

10 The idea was to make a connection to the Type II and Type 
la SNe (not including Type Ib/c), which correspondingly have the 
massive star and compact star origins, respectively. This is however 
not related to the original definitions of Type II and Type I SNs, 
which are based on whether or not there are hydrogen lines in the 
spectrum. 
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pies of Type II and Type I GRBs. We then turn around 
to evaluate the various observational properties (dura- 
tion, hardness, spectral lag, afterglow properties, empir- 
ical correlations, etc) of these Gold Samples, and check 
whether these properties are useful criteria to judge the 
physical category of the bursts. In §6, we discuss the 
intriguing question whether all short/hard GRBs are of 
the Type I origin, and raise the possibility that a fraction 
of (probably even most) high-luminosity short GRBs are 
of the Type II origin. We then dedicate §7 to discuss the 
possible progenitors of GRB 080913 and GRB 090423, 
and argue that most likely they are both of the Type II 
origin. We acknowledge the difficulties of discerning the 
physical origins of GRBs in §8, and cautiously propose 
an operational procedure to associate a GRB with a spe- 
cific model based on multiple observational criteria. Our 
results are summarized in §9. 

2. GRB 080913 AND GRB 090423: INTRINSICALLY 
SHORT/HARD GRBS AT HIGH-Z 

The light curve of GRB 080913 as detected by 
Swift/BAT is shown as the black solid curve in FigQJt. 
The burst duration Tqq (the time interval during which 
90% of the fluence is measured) in the BAT (15-150 keV) 
band is 8 ± 1 s. The average BAT band spectrum can 
be adequately fit by a power law with expo nential cutoff, 
with t he peak energy E p = 93 ± 56 keV (|Greiner et al.l 
I2009al) . A combined Swift/BAT and Konus/Wind (20- 
1300 keV) fit using the Band-function spectrum gives 
E p = 121+§f keV ([Palshin et aT1l2008[) Given the mea- 
sured redshift z = 6.7 (|Greiner et al.l l2009al) . this is 
translated to a rest frame duration of T^ st ~ 1 s, and 
a best-fit rest frame peak energy E r v est ~ 710 keV and 
firest ^ ggg f Qr ^ e cu ^ ff power law and Band- 
function spectra, respectively. Although being recog- 
nized as a long duration burst phenomenologically, this 
burst has an intrinsically short duration and an intrinsi- 
cally hard spectrum. 

In order to compare this burst with other phenomeno- 
logically classified short hard GRBs, we simulate a 
"pseudo" GRB by placing GRB 080913 at z = 1. We 
consider three factors. First, the specific photon flux 
N(E P ) at E p is proportional to (1 + z) 2 /D|, where D L 
is the luminosity distance. This can be translated to an 
increase of a factor of ~ 6.8 of N(E p ) from z = 6.7 to 
2 = 1. Second, we consider the BAT band (15-150 keV) 
emission of the pseudo GRB, which corresponds to an 
energy band lower by a factor of (1 + 6.7)/(l + 1) ~ 3.85 
in GRB 080913. We therefore extrapolate the observed 
BAT spectrum to lower energies and assume a similar 
light curve in that band. Third, we compress the time 
scale by a factor of ~ 3.85 to account for the cosmologi- 
cal time dilation effect. After applying these transforma- 
tions, we are able to construct the BAT-band light curve 
of the pseudo GRB at z = 1 as shown in FigQJi. 

GRB 080913 displ ays a series of early X-ray flares 
(|Greiner et al.ll2009af ). It is interesting to check whether 
they would show up in the BAT band for the pseudo 
GRB to mimic the "extended emissio n" seen in a sub- 
group of Swift "short/ hard" GRBs (jNorris & Bonnelll 
I2006t iTroia et al.1120081, 11 . We therefore manipulate the 
XRT ([Burrows et al 2005b) data of GRB 080913 to sim- 
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Rigorously based on the Tgo criterion, the fraction of Swift 
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Fig. 1. — The simulated 15-150 keV light curves of the pseudo 
GRBs obtained by placing GRB 080913 and GRB 090423 at z = 1. 
The red curves display the extrapolated BAT data, and the blue 
data points show the extrapolated XRT data. Inset: a comparison 
of the light curve of the pseudo GRBs (red) and the observed GRBs 
(black), (a) GRB 080913; (b) GRB 090423. 

ulate the BAT band extended emission of the pseudo 
burst. We first extrapolate the GRB 080913 XRT data 
to the BAT band according to the measured XRT pho- 
ton spectral index. We then follow the three steps men- 
tioned above to shift this BAT-band "virtual" emission 
to the BAT band emission of the pseudo burst. This 
is shown as blue data points in FigOJi. By adding 
the appropriate noise level for the BAT observation, 
we show that these extrapolated XRT emission compo- 
nents stick out the background, which would appear as 
the extended emission in the BAT band for the pseudo 
burst. We note that our method is based on the as- 
sumption of the power law extension of the X-ray flare 

bursts that have Tgo < 2 s is much smaller than that of BATSE 
bursts. Many display extended emission that extends Tgg up to 
several 10s to even more than 100 seconds. The current approach 
in the community is to define a burst "short/hard" if it appears 
short in the BATSE band. A growing trend is to also include some 
bursts with extended emission even in the BATSE band to the 
"short /hard" category. 
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spectrum (0.3 — 10 keV) to the BAT band of the pseudo 
burst (1.3 — 39 keV). On the other hand, since X-ray 
flares are generally b elieved to be due to GRB late cen- 
tral en gine activities dB arrows et alj|20~05at IZhang et al.1 
120061 iLazzati fc Pernal l2007t IChincarini et all l2007f h 
they may have a Band- funct ion or cutoff power law 
spectrum (|Falcone et al.ll2007l ). If the E p 's of the X- 
ray flares are within or not far above the XRT win- 
dow, the extrapolated extended emission would be de- 
graded. We should therefore regard the level of the 
extended emission of the pseudo burst as an upper 
limit. We estimate the BAT-band duration of the pseudo 
GRB as Tgo (pseudo) ~ 2.0 s without extended emis- 
sion or Tgo (pseudo, EE) ~ 140 s with extended emis- 
sion. In any case, the observational properties of this 
pseudo burst are very similar to some "short/hard" 
GRBs detected in the Swift era. By comparing the 
flux level of the pseudo GRB with other short/hard 
GRBs, we find that it belongs to the bright end of 
the short /hard GRB flux distr i bution , similar to, e.g. 
GRB 051221A ([Burrows et all 120061: ISoderberg et al.l 
f2006h. GRB 060313 jlRoming et alJ 12006ft . GRB 060121 
(|Donaghv et~al|[200l . and the recent GRB 090510 de- 
tecte d by Fermi LAT/GBM and Swift (iHoyersten et al 
20091: lOhno et all l2009t iGuiriec et all l2009t iRau et al 



2009) 



GRB 090423 at z = 8.3 is amazingly similar to GRB 
080913. It was detected by Swift / BAT with a BAT-band 
T 90 ~ 10.3 s (jTanvir et all 12009ft. Given t he measured 
redshift z = 8.26±g;g| (|Tanvir et al.ll2009t) . the corre- 
sponding rest-frame duration is ~ Tgo/(l + z) ~ 1.1 s. 
The peak energy measured by BAT is E p = 48.6 ± 6.2 
keV, corresponding to a rest-frame value £J es * = 451 ±58 
keV. We performed a similar analysis on GRB 080913. 
The results are shown in FigJTjD. Nearly identical con- 
clusions can be drawn from both bursts. 

In the above analyses, the intrinsic duration of a 
burst is defined as Tgo/(l + z), and the duration of 
the corresponding pseuodo GRB at z = 1 is defined 
as 2Tgo/(l + z). These calculated durations correspond 
to different energy bands in the rest frame (the same 
observed band after redshifting) . Strictly speaking, in 
order to derive the durations of the pseudo GRBs in 
the observed energy band, one needs to know the time- 
dependent spectral information, which is not available 
for these bursts. Observationally, pulse widths at high 
energ ies tend to be n arrower than those at low ener- 
gies jFord et all Il995t iRomano et aTl 120061: iPage et all 
l2007f) . An empi rical relation w oc E~ a with a ~ 0.3 has 
been suggested dFenimore et al.lll995t iNorris et al.ll2~005l : 
iLiang et all 12006ft . For' a given observed energy band, 
this suggests woe (l + z) _a , which would correspond to 
a correction factor of (1 + z) a_1 rather than (1 + z) -1 
to derive the intrinsic duration. However, GRB prompt 
emission is usually composed of multiple pulses. The 
separations between the pulses, which are more relevant 
for the Tgo definition, may not follow the same energy- 
dependence of the pulse widths. We therefore do not 
introduce this extra correction factor of Tgo throughout 
the paper. For GRB 080913 and GRB 090423, if one 
takes the (1 + z) a ~ Y correction factor, the derived in- 
trinsic durations are in the marginal regime between the 
phenomenologically-defined long and short GRBs. 
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Fig. 2.— The T go - HR diagram of GRBs. The background 
orange dots are BATSE GRBs. Overplotted are Type II Gold 
Sample (blue), Type I Gold Sample (red), and other short/hard 
GRBs (green), mostly detected by Swift. Open symbols are for 
the observed values, while the filled symbols are the rest-frame 
values. For short GRBs with extended emission, those with the 
short spike only are denoted as circles, while those including the 
extended emission are denoted as squares. The same bursts (with 
different Tgo with or without extended emission) are connected by 
lines. GRB 080913, GRB 090423, their pseudo counterparts at 
2 = 1, and their rest-frame counterparts are marked with special 
colors / symbols. 

Figure [3 displays the locations of GRB 080913, GRB 
090423, their corresponding pseudo GRBs at z = 1, 
and their rest-frame counterparts in the traditional 
Tgrj— HR (hardness ratio) two-dimensional distribution 
plane. Also plotted are the BATSE GRB sample (or- 
ange), the Gold samples of Type II (blue) and Type I 
(red) GRBs, and the Other SGRB Sample (green) (see 
§5 for the details of the sample definitions). It is evident 
that GRB 080913 and GRB 090423 would have been rec- 
ognized as phenomenologically short/hard GRBs should 
they have occurred at z < 1. 

3. PHENOMENOLOGICAL VS. PHYSICAL 
CLASSIFICATION SCHEMES: WEAKNESSES AND 
STRENGTHS 

The eventual goal of GRB studies is to identify the 
physical origins of every observed GRB, including its pro- 
genitor system, central engine, energy dissipation mech- 
anism, and radiation mechanism. To achieve this goal, 
a combination of observations and thcorcical modeling is 
needed. The number of competive models and the al- 
lowed parameter space steadily reduce as more and more 
observational data are accumulated. This is evident in 
the history of GRB studies: while more than 100 mod- 
els were proposed before 1992 (|Nemiroff 1 11 994). only two 
broad categories of progenitor models remain competi- 
tive at the time when this paper is written. A group of 
GRBs are hosted by active star-forming dwarf galaxies 
(jFruchter et al.lr2006l ). s ome of which have clear (Type Ic) 
super nova associations dHiorth et ah 2003; St anek et al.1 
120031: ICampana et"aH 120061: iPian et all 12006ft ! This 
points towards a massive star origin of this group of 
bursts. At least a few bursts were discovered to be as- 
sociated with galaxies with a very low star forming rate 
(jGehrels et al.l l2005: Bl oom et al]|2006t iBarthelmv et all 
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l2005al iBerger et al.l I2005D , which point towards a non- 
massive-star origin of the bursts, likely due to mergers of 
compact objects. Therefore it is now justified to discuss 
at least two physically distinct categories of GRB models 
as well as how to associate a particular burst with cither 
category based on certain observational criteria. 

In the literature, some physical classification s chemes 
of GRBs have b een discussed (jZhang et al.l l2007bl : 
iBloom e~a l. 2008). Strictly speaking, these are not clas- 
sifications of GRBs, but are classifications of models that 
interpret GRB data. A scientific classification scheme is 
based on statistical formalisms, which make use of a uni- 
form set of observational data with instrumental biases 
properly corrected, and classify objects based on statisti- 
cally significant clustering of some measured properties. 
Examples include to classify supernovae broadly into 
Type H/I based on whether there are/are not hydrogen 
lines in the optical spectrum , and to cl assify GRBs into 
two ( iKouvelio tou et al.lfl993l ) or three (jMukhcric e et al.l 
I1998L iHorvath Ifl998h classes based on BATSE T 90 anal- 
yses. The classes defined by the phenomenological data 
do not carry physical meanings, and theoretical model- 
ing is needed to clarify whether different phenomenolog- 
ical classes of objects are of different physical origins. 
Compared with the SN classification schemes, which are 
based on the "yes/no" criteria regarding the existence 
of spectral lines and therefore are relatively insensitive 
to the instrumental details, the GRB phenomenolog- 
ical classification schemes suffer another major draw- 
back, i.e. every parameter that one can directly mea- 
sure is strongly instrument-dependent. For example, T90 
is strongly energy-dependent, and sensitivity-dependent, 
so that a "short" GRB in a hard energy band would be- 
come a "long" GRB in softer bands or if the detector 
sensitivity is increased. The membership of a particular 
GRB to a particular category (e.g. long vs. short) is 
not guaranteed. As a result, such classification schemes 
cannot be compared from one mission to another, and 
are of limited scientific value. 

A physical classification scheme, on the other hand, 
is on theoretical models that interpret the data. As a 
result, it suffers the great difficulty of associating a par- 
ticular burst to a particular model category. In order to 
achieve the goal, multiple observational criteria are de- 
manded, but always with non-uniform instrumental se- 
lection effects. Ideally, with infinitely sensitive detectors 
in all wavelengths, it may be possible to derive a set 
of quantitative observational criteria that can be used 
to rigorously associate a particular GRB to a particular 
model category based on statistical properties. However, 
realistically this is essentially impossible since different 
criteria rely on completely different observational instru- 
ments with different observational bands and sensitivities 
which are quite non-uniform. Also different criteria could 
carry different weights in judging the associated model 
category of a particular burst. The weighting factors of 
different criteria are also difficult to quantify. Human 
insights rather than pure statistical analyses are needed. 
Another drawback of a physical classification scheme is 
that it depends on the models, which are subject to fur- 
ther development as more data are accumulated. The 
classification criteria are therefore also subject to modi- 
fication based on data. This can be diminished by invok- 
ing model-independent criteria as much as possible. For 



example, the Type I /II GRB model classification scheme 
discussed in this paper only appeals to whether the model 
invokes a degenerate-star or a massive-star, regardless of 
the concrete progenitor systems or energy dissipation and 
radiation mechanisms (see §4 for full discussion). 

Despite of its weaknesses, a physical classification 
scheme of models and associating a particular object to 
a particular model class has the strength to achieve a 
better understanding of the physical origin of astrophys- 
ical objects. For example, in the supernova field, there 
is now a consensus that only a sub-group of Type I SNe 
(Type la) has a distinct physical origin, which is related 
to explosive disruptions of white dwarfs. The other two 
sub- types of Type I SNe (Type Ib/Ic) are more closely 
related to Type II SNe and form together a broad phys- 
ical category of SN models that invoke massive star core 
collapses. Such a physical classification scheme of SN 
models (massive star origin vs. white dwarf origin) and 
the efforts to associate the observed SNe to them reflect 
a deeper understanding of the physical origins of SNe. 
The same applies to GRBs. The statistical classifica- 
tion of long-, short- and probably intermediate-duration 
GRBs has been established in the BATSE era. However, 
it took several missions and many years of broad-band 
observations to reveal that there are at least two phys- 
ically distinct types of models that are associated with 
these GRBs. Although data are not abundant enough to 
unambiguously associate every individual GRB to these 
model categories, current data already revealed some 
perplexing observational facts (§1) that demand more se- 
rious investigations of the observational criteria to judge 
the physical origin of a particular GRB (i.e. the physical 
model associated with this GRB). 

In the rest of the paper, we will discuss Type I/II 
GRBs, which are defined as the GRBs that are asso- 
ciated with two distinct physical models. This is not a 
new classification scheme of GRBs to replace the exist- 
ing long/soft vs. short/hard classification scheme, but is 
a parallel classification of the models that the observed 
GRBs can be associated with based on multiple crite- 
ria data analyses. The two approaches are complemen- 
tary. As discussed above, Tqq is energy-band-dependent 
and sensitivity-dependent, so that the membership of a 
particular GRB to a particular duration category is not 
always guaranteed. On the other hand, if adequate infor- 
mation is retrieved in an ideal observational campaign, 
the association membership of a particular GRB to a 
particular physical model category is almost certain re- 
gardless of the detector energy band and sensitivity. For 
example, if a SN is detected to be associated with a GRB, 
one can safely associate this GRB to the Type II model 
category regardless of its Tg detected by different detec- 
tors. 

4. TYPE I/II GRBS, THEIR OBSERVATIONAL CRITERIA, 
AND PHYSICS BEHIND 

We reiterate here the definitions of the Type I/II 
GRBs. Imp r oving upon the descriptions presented in 
I Zhang et al.l (|2007b[ L we hereby more rigorously define 
the following: 

• Type I GRBs (or compact star GRBs) are those 
GRBs that are associated with the theoretical 
models invoking destructive explosions in old- 
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population, degenerate, compact stars. The like- 
liest model candidate is mergers of two compact 
stars. 

• Type II GRBs (or massive star GRBs) are those 
GRBs that are associated with the theoretical 
models invoking destructive explosions in young- 
population massive stars. The likeliest model can- 
didate is core collapses of massive stars. 

Here we do not specify the progenitor systems of 
each model type. In reality, there could be mul- 
tiple possible prog enitor systems wit hin each model 
category (see also iBloom et al.l 120081 ). Within the 
Type I model category, possible progenit o r systems in- 
clude NS-NS mergers (IPaczvnskH 119861: lEichler et all 



119891 : INaravan et "a l. 1992; Ross wog et aTlj2003h . NS-BH 
mergers (jPaczynskil 119911 : iFaber et al.l |2006|L and pos 



sibly BH-WD o r NS -WD mergers ( Fryer et alj 119991: 
iKing et all l2007l)(c.f. INaravan et alj hoOll ). see iNakarl 
(|2007| ): lLee fc Ramirez-Ruid (|2007| ) for reviews. On the 
other hand, within the Type II model category, one may 
have collapses of single stars (i.e. collapsars, Woosley 
1993; MacFadyen & Woosle y 1999), or collaps es of mas- 
sive stars in binary systems (|Frver et al.ll2007| ). 

The definitions of Type I /II GRBs are based on the 
physical models that GRBs can be associated with rather 
than their observational properties. The scheme is there- 
fore intended to be "operational" . The connections be- 
tween the physical model properties and the observa- 
tional criteria are not straightforward, and probably very 
difficult for some GRBs. 

In the following, we review a list of observational cri- 
teria discussed in the literature that may be applied to 
differentiate the two physically distinct model categories 
that GRBs can be associated with (e. g. iDonaghv et al.l 
120061: IZhang et al.ll2007bt IZhangj|2006[ ). and discuss the 
physical justifications of each criterion. As justified be- 
low, some criteria (e.g. §3.1-§3.3) are more directly 
related to the progenitor system of a GRB. On the 
other hand, the traditional criteria invoking the observed 
gamma-ray properties (e.g. §3.4-§3.7) are more related 
to radiation physics and have less direct connection with 
the progenitor system. Some afterglow properties (§3.8, 
§3.9) do carry information of the progenitor, but theoret- 
ical modeling is invoked (and, hence, less definitive be- 
cause of the uncertainties inherited in the models). The 
statistical properties (§3.10, §3.11) can be related to the 
progenitor system, but again large uncertainties are in- 
volved in the identification of the explicit progenitor sys- 
tem and its cosmological evolutionary scenario. The best 
clue may be gravitational wave signals (§3.12). However, 
they are beyond the current detector capability. 

4.1. Supernova Association 

A positive detection of a supernova (SN) signature as- 
sociated with a GRB would undoubtedly establish the 
association of the burst with Type II. However, the sam- 
ple of the robust GRB-SN associations is currently small. 
Non-detections of a SN signature could be due to multi- 
ple reasons, e.g., the afterglow is too bright so that the 
SN light is buried beneath the afterglow level; the fol- 
low up observations were not "deep" enough or not at 
the right time window; or the lack of an underlying SN 



is genuine. Only the last case is helpful to judge the 
physical model category of a burst, although the con- 
clusion is still not clear cut. A genuine SN-less GRB is 
certainly consistent with the Type-I origin. However, it 
has been discussed in the literature that some core col- 
lapse GRBs m a y not eject enough 56 Ni to power a SN 
(IWooslevi 119931: [Heger et al.l 120031 iNagataki et "all 120031 . 
120061: iTominaga et al.ll200'^ v! so that the lack of a genuine 
SN signature may not be evidence completely against the 
Type II origin. On the other hand, we notice that the 
large uncertainties involved in SN explosion physics pre- 
vent the models from having a definite predictive power 
regarding the SN signature. Looking back into the his- 
tory, the predictions of the SN signature accompanying 
GRBs have followed a serpentine (and ironic) path. The 
first co re-collaps e GRB model was dubbed "failed super- 
nova" (|Wooslevlll993h . which predicts no SN signature 
associated with a GRB. Driven by the possible GRB 
980425/ SN 1998bw association, the model was developed 
to allow a SN associated with a GRB within the "col- 
lapsar " scheme. According to MacFad ven fc Wooslevl 
(1999), the model predicts that "collapsars will always 
make supernovae similar to SN 1998bw". Indeed the 
statement that "the data and models are consistent with, 
though not conclusive proof of, the hypothesis that ALL 
long-soft GRBs are accompanied by SNe of Type Ic" was 
made ri ght before the discovery o f GRB 060614 and GRB 
060505 (|Wooslev fe Bloorn1l2006D . Would the discovery of 
the SN-less long GRBs (060614 and 060505) then beg for 
a dichotomy of core-collapse GRBs (one group with and 
another group without the SN association)? Although 
this is certainly plausible, a simpler picture would be 
that all genuine SN-less GRBs have the Type I origin. 
In this paper, we take lacking a genuine SN as a support 
to the Type I GRB, but do not take this criterion alone 
to define a Type I GRB. On the other hand, since there 
is no observational fact that demands the existence of 
SN-less Type II GRBs 12 , we do not automatically asso- 
ciate any genuine SN-less long GRB with the Type II 
(or Type II candidate) physical model categories unless 
there are other strong supports to the scenario (see §4.1 
and §7 for details). 



4.2. Star Forming Rate of Host Galaxy 

Type II GRBs are related to massive star deaths, so 
they must reside in host galaxies with active star for- 
mation. So star forming rate (SFR), or more rigorously, 
specific star forming rate (SSFR, i.e. SFR per unit mass) 
of the host galaxy is a critical parameter to judge the 
membership of Type II GRBs. 

On the other hand, compact star mergers can occur in 
host galaxies both with and without active star format ion 
(jBelczvnski et al.ll2006t IZheng fc Ramirez-Ruizll2007l) . If 
we see a GRB residing in an elliptical or an early type 
galaxy, we are more confident that no massive star is 
involved in the event, and that the burst should be as- 
sociated with Type I. Those Type I GRBs residing in 
star-forming galaxies are more difficult to identify. Due 

12 In our opinion, GRB 060614 and GRB 060505 are not 
solid Type II cand i dates. As will be discussed in §5 (see also 
IGehrels et a l. 2006; Z hang et al1l2007bl) . GRB 060614 is likely a 
Type I GRB. GRB 060505 has a much lower overall energetics (in- 
cluding gamma-ray and afterglow) than most other Type II GRBs. 
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to the additional time delay required for the two com- 
pact objects to coalesce, a Type I GRB site is expected 
to be more aged than the site with active star formation. 
As a result, at least some Type I GRBs should preferen- 
tially reside in the regions with relatively low SSFR in 
the star forming host galax y. On the other hand, there 
are channels of fast mergers (jBelczvnski et al.1l2006l ) that 
lead to almost "prompt" mergers of compact stars. In 
such a case, Type I GRBs can reside in high SFR regions 
within star-forming galaxies. 

4.3. Position Offset with Respect to Host Galaxy 

A related criterion is the offset of the GRB location 
with respect to the center of the host galaxy. The phys- 
ical motivation is that Type I GRBs invoke mergers of 
binaries including at least one NS, which likely received 
a "kick" at birth so that the binary system would mi- 
grate from its original birth location. By the time when 
the two compact stars coalesce, the system should have 
a large offset from the galaxy center or even be out- 
side of the host galaxy (Bl oom et al.lll999h . Indeed sev- 
eral Gold Sample Type I GRBs show such a property 
dGehrels et al J 12005 Bloom et alJl2006t iFox et al.ll2005l; 
Berger et all 120051: fearthelmv et all l2005at iTroia et all 
20081 ). On the other hand, Type II GRBs explode right 
at the location where the progenitor stars are formed 
and therefore should b e in the star formi ng regions in- 
side the host galaxy (|Bloom eTal1l200l . This is m 
general consistent wi th the observations of long GRBs 
(|Fruchter et al.l [20061 ). Outliers do exist. For example, 
GRB 070125 is a long GRB whos e birth location is in a 
galactic halo ijCenko et al.ll2008bf ). 

Complications arise if a GRB is found not inside any 
galaxy. It is difficult to judge whether a GRB is "kicked" 
out from a nearby host galaxy whose projected image is 
near the location of the GRB, or it is associated with 
a more distant galaxy at high-z. This problem arises 
for a good fraction of short/hard GRBs detected in the 
Sw ift era. For e xamp le, GRB 060502B was suggested 
by IBloom et all (|2007h to be associated with a nearby 
galaxy at z = .287 (with a l arge a ngular offset), while 
it is included by [Berger etai] (2007) as one of the high-z 
missing- host short/hard GRBs. 

4.4. Duration 

Theoretically, we do not know exactly which time scale 
defines the GRB duration. In principle there are three 
time scales that are relevant. The first one is the duration 
of the central engine activity t en gine . This corresponds to 
the accretion time scale of an accretion-powered central 
engine model (usually invoking a black hole - torus sys- 
tem) , or the spindown time scale of a spindown-powered 
central engine model (usually invoking a rapidly rotat- 
ing millisecond magnetar or a maximally rotating black 
hole whose spin energy is tapped via a magnetic torque 
through the Blandford-Znajek mechanism). The second 
time scale is the time scale tj e t during which a relativistic 
jet is launched. In principle there could be epochs dur- 
ing which a jet is launched, but it is not relativistic or 
not relativistic enough to power the observed gamma-ray 
emission. The third time scale is energy dissipation time 
scale tdis- Current Swift observatio ns suggest that th e 
GRB prompt emission is "internal" ( Zhang et al.l 12006) . 
This requires that the relativistic jet dissipates energy 



internally before being decelerated by the external cir- 
cumbu rst medium. The dissipa tion could be via internal 
shocks (|Rees fe Meszarosil994l) or magnetic reconnection 
(|Usovlll992[ iThomp son 1994). In principle, one can have 
an active central engine without launching a relativistic 
jet, or have a relativistic jet without significant dissipa- 
tion. In general, the observed GRB duration Tgo (which 
also depends on the energy band and the sensitivity limit 
of the detector) should satisfy 13 

Tgo 5^ tdis tjet — tengine ■ (1) 

In most studies, however, T 90 ~ t en gine ~ t jet ~ t dls has 
been assumed. 

If Tgo is equal to or at least is proportional to t en gine , 
as is assumed by most central engine modelers, then the 
duration information may be tied to the progenitor prop- 
erties of GRBs. In particular, Type II GRB progeni- 
tors have a massive envelope, which can power a long- 
duration GRB through accretion. According to the col- 
lapsar scenario (|MacFadyen fe Wooslevl[l999 ). the dura- 
tion of the burst is defined by the envelope fallback time 
scale, which is typically 10s of seconds. The model there- 
fore suggests that Type II GRBs should typically have 
long durations. On the other hand, NS-NS and NS-BH 
m ergers typically h ave an accretion time scale ~ 0.01-0.1 
s (jAloy et al.l 12005) if the central engine is a BH-torus 
system. Therefore Type I GRBs should typically have 
short durations. Indeed, a 1-second duration burst is al- 
ready too long to be accommodated within the simple 
merger scenarios. One needs to introduce additional in- 
gredients (e.g. an interm ediate neutron star p hase) to 
increase the duration (e.g. IRosswog et al]|2003f ). 

This physically-motivated clear dichotomy was broken 
in the Swift era. Swift discovered that X-ray flares pre- 
vail in more th an half GRBs, in both long and short 
duration ones (jBurrows et all l2005at IChincarini et all 
120071: iFalcone et alll2007f ). This suggests that the GRB 
central engine activity is not limited to the prompt 
phase, and is m uch longer than T 9 q in both long 
and short GRB s llZhang et all [20061 iFan fc Weil [2001 
llazzati fc Pernall2007D . The progenitor and central en- 
gine models must then be modified to invoke a much 
longe r accretion time scale llKing et al . 2005; Pcrn a et all 
2006; Pro gafc Zhanel [27X1611 or a non-BH-torus central 
engine (|Dai et al.l 120061: IStaff et all l2007f ). More im- 
portantly, several strong Type I GRB candidates (e.g. 
GRB 050 724) are not short, but have softer, extended 
emission (IVillasenor et alll2005l: iBarthelmv et al.H 2005a: 
iNorris fc BonneU I2006T) . " The merger models therefore 
must be modified t o account for this extended emis- 
sion (|Rosswogll2007l ). Type I GRBs no longer must be 
"short" . 

The discussion above only applies to the case when 
the line of sight pierces into the relativistic jet, i.e. 
the on-beam geometry. In this case, the observed 
time scales reflect the time scales at the central engine 
(Kobavas hi et al.lll997l ). In the case of an off-beam geom- 
etry, i.e. the jet with opening angle 9j is beaming towards 
an angle 6 V > 0j with respect to the line of sight, the 
observed time scale no longer traces that of the central 

13 Here we have assumed that Tgo records the GRB internal 
emission only. This is true for most cases. Occasionally the ob- 
served prompt emission may also include the emission from the ex- 
ternal shocks. T90 should be removed from Eq.JTJ for these cases. 
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engine. For a discrete pulse, if the pulse duration solely 
reflects the duration of the emission powered by the cen- 
tral engine, i.e. the rising and falling of the lightcurve 
reflects the increase and decrease of the central engine 
luminosity (in constrast to those models that interpret 
the decaying wing as the high-latitude emission) , the ob- 
served duration off beam is related to the on-beam value 
through the ratio of the Doppler factor (given the same 
comoving value), i.e. 

t(oS beam) _ V(0 = 0) 1-/3 cos(6»„ - 9j) 

i(on beam) ~ V{9 = 9 V - 9 3 ) ~ 1-/3 ' 

(2) 

where the Doppler factor is defined by 

_ 1 

v = r(i-/?cos0) ' (3) 

and 9 is the angle between the line of sight and the veloc- 
ity vector of the ejecta, which is taken as the closest ap- 
proach to the jet (9 v — 9j). For multiple emission episodes 
(i.e. multiple pulses in the light curve), the time interval 
of the quiescent episodes do not vary with the viewing 
direction. So Eq.@ applies to the total duration of a 
GRB only if the prompt emission has one single pulse. 
Also since the observed flux is lower for a lower T> 1 given 
a same detector sensitivity, the off-beam Tqq tends to be 
shorter than that predicted by Eq.(2) due to the limiting 
flux threshold effect. 

The off-beam model predicts that the afterglow light 
curve should display a rising beh avior initially before the 
1/r beam enters the line o f sight ( Panaitc scu fc Meszaros! 
119991 : iGranot et al.l l2002| ). Broadband observations of 
the majority of GRB afterglows do not show such a sig- 
nature. So the off-beam geometry, if any, is rare. 

4.5. Hardness 

The connection between the hardness of spectrum and 
the GRB progenitor is less direct. It is related to the 
unknown internal energy dissipation mechanism and ra- 
diation mechanism, which in turn depends also on the 
composition of the GRB ejecta. GRB spectra are usu- 
ally categorized as a smoothly-joined power-law, namely, 
the Band- function (Ba nd et al.l Il993 h . The hardness of 
a GRB is likely related to the location of the E p , but 
the flatness of the spectral index below E p may also 
play a role. Theoretically, the spectral slope is more 
close ly related to the particle acceleration mechanism 
(e.g. ISironi. fc Spitkovskv II2009T) and the "co mpactness" 
of the emission region (e.g. lPe'er et a Tl l2006D . The spec- 
tral peak energy, E p , can be related to the GRB emis- 
sion model parameters more directly , although model- 
dependent ([Zhang fe Meszaros 20023). We will mainly 
discuss the E p models more closely in the following. 

In general, E p is a function of the burst luminosity L, 
the Lorentz factor T of the ejecta, and the radius R of the 
emission site from the central engine. In order to address 
whether a GRB is hard or soft, one needs to specify a 
particular emission model. In the following we discuss 
three internal emission models currently discussed in the 
literature. 

Internal shock model. Within this model, the 
gamma-ray E p can be defined either by synchrotron ra- 
diation or synchrotron self-Compton (SSC). In general 
one can write E p ~ Thj^eB' /mc), where k — (2,4) 



for synchrotron and SSC, respectively. Since the comov- 
ing magnetic field strength in the ejecta flow satisfies 
B' oc L 1 / 2 i? _1 r _1 for both t he ordered and the random 
magnetic field components (jZhang fc Meszaroa l2002al ) . 
one has 

Ef cx 7e fc L 1 / 2 iT 1 (l + z)- 1 oc 7e fc J L 1 / 2 r- 2 ( 5^ 1 (l + z)- 1 , 

where L is the initial kinetic luminosity of the ejecta, 
St is the variability time scale of the unsteady GRB 
ejecta wind, and the internal shock radius is R ~ T 2 St. 
Note that E p is negatively correlated with the bulk 
Lorentz factor (oc T -2 ), which is contrary to the intu- 
ition that high T bursts should be hard. Here 7 e is the 
characteristic Lorentz factor of the electrons that con- 
tribute to the emission at E p . Under the fast cool- 
ing condition, which is generally satisfied for internal 
shocks, 7 e corresponds to the minimum "injection" en- 
ergy of the electrons, which is related to the "relative" 
Lorentz factor between the two colliding shells Tf s , i.e. 
7 e oc r /s ~ (Tf/T s + T s /T f )/2, where T f and T s are the 
Lorentz factors of the fast and slow shells, respectively. 
If the r variation of a flow is proportional to the average 
Lorentz factor T, i.e. Ar oc T or Tf/T s ~ const, then 
7e essentially does not depend on F, so that a higher E p 
should correspond to a lower T. On the other hand, it 
is possible_that high-r flows may be more variable, e.g. 
F//F s oc T. If this is the case, then the negative de- 
pendence on r in Eq.((4]) is canceled out (for k — 2) or 
reversed (for k = 4). In the traditional internal shock 
model, the variability time scale St of the ejecta can be 
derived from the observation. Analyses of power density 
spectra of GRB light curves (jBeloborodov et al"1ll998f) 
suggest that the GRB temporal behavior may be self- 
similar, i.e. lacking a characteristic time scale. In the 
past, the minimum variability time scale, which can be 
as small as milliseconds for both short and some long du- 
ration GRBs, has been adopted to estimate the internal 
shock radius. Alternatively, it is possible that the rapid 
variability in GRB light curves may be caused by other 
mechanisms, su ch as relativistic turbule nce inside the 
emission region (jNaravan fc Kumar! [20091 ) . Within this 
latter scenario, the outflow variability time scale relevant 
to internal shocks can be much longer. Physically, Type I 
GRB outflows may directly carry the variability informa- 
tion from the inner central engine, i.e. the dynamic time 
scale of the innermost accretion torus around the black 
hole, St ~ t dyn ~ l2V3ir{GM bh /c 3 ) ~ l(M bh /3M Q ) 
ms (where M^h is the mass of the black hole), or the 
spin period of the central magnetar or black hole, St ~ 
Pengine ~ 1 m S- On the other hand, a Type II jet needs 
to pentrate through the heavy stellar envelope so that 
the initial temporal information from the inner central 
engine may be smeared out and regulated. The observed 
variability time scale may be related to that of fluid in- 
stabilities, and therefore could be much longer. If abun- 
dant pairs are produced, it has been argued that the pair 
photosphere would effectively screen o ut the variability 
time scales smaller than a c ritical value (jKobavash i et al.l 
2002; Meszaros ct al. 2002). This is more relevant to low- 
r events for which the internal shock radii are below the 
pair photosphere. 

With all these complications in mind, one may compare 
the expected E p for Type I and Type II GRBs based 
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on Eq.(|4]). On average, Type I GRBs have an isotropic 
gamma-ray luminosity L 2-3 orders of magnitude smaller 
than that of Type II GRBs (see the theoretical argument 
in §3.8, and the observational data in §4 and Table 1 
below). On the other hand, St of Type I may be smaller 
than that of Type II by 2-3 orders of magnitude. This 
gives 



4 s '(I) 
MS(ii) 



(10-30) 



[ 7e fc r- 2 (i + *)-!](!) 
[ 7e fc r-2(i + *)-i] (ii) 



(5) 



This suggests that in a large parameter space Type I 
GRBs can be harder than Type II GRBs. If 7 e is similar 
for both types, Type I GRBs can be harder than Type II 
GRBs as long as their bulk Lorentz factors are not larger 
than those of Type II by a factor more than (3 — 5) times. 
If 7e oc r, Type I are generally harder than Type II for 
both the synchrotron model (k — 2, regardless of the 
value of r), and the SSC model (fc — 4, the E p ratio is 
positively dependent on T. Theoretically, Type I GRBs 
should have higher T's due to their less baryon loading 
as compared with Type II GRBs. This favors a harder 
spectrum of Type I even more for the SSC model. A 
systematically smaller redshift z for Type I GRBs (due 
to the merger delay with respect to star formation) also 
helps to increase the hardness contrast between the two 
types. In reality, there are large dispersions in L, St, 
r, 7 e and z in both types. On the other hand, the HR 
distribution of the BATSE short/hard vs. long/soft di- 
chotomy also shows a large dispersion (Fig[2|). In general, 
the statement that Type I GRBs are harder than Type 
II GRBs can be made within the internal shock models 
in the statistical sense. For individual bursts, one can- 
not draw a firm conclusion regarding the hardness of a 
particular burst due to the large uncertainties involved 
in the parameters. 

Photosphere model. The possibility that the 
observed GRB emission has a do minant contribu- 
tion from the fireball photosphere ([Thompson! |1994t 
Me szar os fc Reesll2000t IMeszaros e t al. 2002) has gained 
increa sing attention rec ently (iRees fc MeszaroL 
Rvdd 120051: iPe'er et all 120071: iThompson et al.l 



2005 



2007 



Ioka etalJ 120071: IGhisellini et all 120071: iRvde fc Pe'erl 
20081 lLazzati et al.l l2009h . In this model, E p is related 
to the observed photosphere temperature T p h- For a 
"naked" fireball, i.e., a fireball expanding into a vac- 
uum, the observed photosphere temperature (and hence 
E p ) depends on whether the photosphere radius R p h is 
below or above the fireball coasting radius R c . For a 
large dimensionless entropy of the fireball rj > r\ c i ~ 
10 4 [i52^?o y] 1 ^ 3 ( wri ere i?o is the initial radius of the fire- 
ball. Throughout the text the convention Q n = Q/10 n is 
adopted in cgs units.) 14 , the fireball becomes transparent 
during the acceleration phase (i.e. R p h < R c )- The ob- 
served fireball temperature is essentially the temperature 
at the central engine, i.e. T p h ~ To, so that 

Ef' 1 ~ kT (l + z)- 1 ~ L^R-^il + z)- 1 . (6) 

Thi s is the regime discussed in most photosphere mod - 
els l)Thompsorili99l IMeszaros fc Reedl20"OOtlRvddl2005l) . 



14 This critical entropy is derived (Zhang & Mcszaros 2002a) 
within the discrete shell regime, rather than the con tinuous wind 
regime IjMeszaros & Roes 2000; Mcszaros ct al. 2002). This is usu- 
ally justified, since typically one has rj > r\ c \ in this regime. 



On the other hand, if the fireball becomes transparent 
beyond the coasting radius (R p h > R c ), the photosphere 
temperature drops with radius due to the decrease of 
residual internal energy during the expansion, so that 
T ph = T (R C /R ph ) 2 / 3 (|Meszaros fc Reesl |2000[) . The 
detailed parameter dependences are related to whether 
the opacity is de nned by a discrete sh ell or a continu- 
ous outflow wind (|Meszaros et al.ll2002D . For the former 

(Vd < V < Vc2, where r\ cl ~ 250[L 5 2i?o~7] 1/5 ), one has 
(|Zhang fc Meszarosll2002al ) 



ctM ~L- l ' X2 R- 1/6 T{l + zY 



(7) 



(8) 



For the latter (77 < ?/ c i), one has 

Ef > 3 ~ L- B /"ijj/«r8/3(i + z )-i . 

If additional energy dissipation occurs at small radii, pair 
production can occur which enhanc es photon opacity and 
increases the photosph ere radius (|Meszaros et al.l 120021 : 
iRees fc Meszarosll2005l ). We note that the "naked" fire- 
ball scenario is more relevant to Type I GRBs. 

With the presence of a stellar envelope, the photo- 
sphere emission of a Type II GRB is likely modified. 
Due to continuous energy dissipation and heating inside 
the envelope, the jet cannot reach the maximum Lorentz 
factor but instead stores a significant energy in heat be- 
fore erupting out from the envelope. Effectively, the 
GRB fireball "central engine" is moved from the central 
black hole or ma gnetar to th e loc ation slightly belo w the 
stellar envelope (IThompsonll2006l:lThompson et alj|2007t 
IGhisellini etaf] 12007( 1 . The jet at this radius R ~ i?» 
has a moderate Lorentz factor and a comoving tem- 
perature Tl ~ (L/AnT^R^a) 1 / 4 , and an observer frame 
temperature T* = r»T^. As the jet erupts out from the 
envelope, it will undergo rapid acceleration under its own 
thermal pressure. If i?» is greater than photosphere ra- 
dius for a naked central engine, the fireball would become 
transparent shortly after exiting the star due to the rapid 
fall of density. So the real photosphere radius is essen- 
tially R p h ~ and R p h < R c is always guaranteed. 
The peak energy E p is defined by T p h — T*. This leads 
to a variation of Eq. ([6]) in the form of 



m>h,l 
p 



L^T l J 2 R- l '\l + z)- 



(9) 



Within the photosphere models, it is not obvious why 
Type I GRBs should be systematically harder than Type 
II GRBs. The trend, if any, should be opposite. The 
logic is the following. First, given the same parameters 
of L, R and z, o ne typically has £g M > E ^ 2 > E^ 3 
(e.g. Eq.(24) of IZhang fc Meszarod r2002a( ) . Next, the 
stellar envelope effectively "raises" the photosphere for 
Type II GRBs, so Eqs.© and © are usually not rele- 
vant. One therefore may only compare the case of Eq.([6]) 
for the two types of GRB, since Eq.© can be related to 
Eq.© through i? = R*/T*. Equation © suggests that 
Type I GRBs, typically with smaller L, should be softer 
than Type II GRBs at the same redshift. A smaller z for 
Type I GRBs may compensate their softness, but in gen- 
eral it is not straightforward to claim that Type I GRBs 
should be systematically harder than Type II GRBs for 
the photosphere model. Pairs may lower the photosphere 
temperatures of some high-L GRBs (especially for Type 
II), which may help to account for the observed trend, 
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but no handy analytical formula is available to perform 
direct comparisons. 

The recent Fermi-detected GRB 080916C (|Abdo et al.l 
I2009D showed a series of featureless Band-function spec- 
tra. The expected photosphere emission component is 
missing, suggesting a Poyntin g flux dominated flo w at 
the base of the central engine (|Zhang fe Pe'enl2009h . At 
least for this burst, the observed E p is not the thermal 
peak of the photosphere emission. 

Magnetic dissipation model. Finally, if the 
GRB 
1991 



outflow is Poyn ti ng fl ux dominated 
iMeszaros fc Reesl 119971: ISpruit et al.l 



Lvutikov fc Blandford 2003; Liang & NoguchU l2009f ) 



the characteristic frequency of emission would take a dif- 
ferent form and have different dependences on the ejecta 
parameters. The locations of the magnetic reconnection 
regions are unknown. If dissipation occurs at small radii, 
the effect is to modi fy the photosph ere emission through 
continuous heating (|Gianniosll2008D . This is effectively 
a photosphere model, which has been discussed above. 
Alternatively, a Poynting-flux-dominated outflow can 
reach a global dissipation at a large radius where the 
MHD approximation is broken (Usovlll994t ISpruit et al 



2001 



2003). 



IZhang fc Meszarosl l2002at iLvutikov fc B landford 



Lacking a macroscopic reconnection model, E p in the 
reconnection model is difficult to calculate. Under differ- 
ent assumptions, the expression of E p may take different 
forms. For example, in a random electric/magnetic field 
in the magnetic reconnection region, electron accelera- 
tion may be balanced by radiation cooling. The typical 
electron Lorentz factor is therefore j e <x B~ x l 2 . The 
synchrot ron peak energy may be then expressed in the 
form of (|Zhang fc Meszarodl2002ah 

E^ 9 oc T(l + z)- 1 , (10) 

which depends on T and z only. Type I GRBs can then 
be harder than Type II GRBs, again because Type I 
GRBs tend to have a cleaner environment, and hence, 
less baryon loading, than Type II GRBs. 

Similar to the duration discussion (§4.5), the above 
discussion applies to the on-beam geometry. For an off- 
beam geometry, the observed Ep is smaller by a factor 
of the Doppler fa ctor ratio (Eq.[2] ) . This effect has been 
discussed by, e.g. lYamazaki et al.l ((2004b). 

4.6. Spectral Lag 

Soft GRB emission usually arrives later than hard 
emission in some GRB s. This "spectral lag" is evident for 
long-duration GRBs (iNorris et alJ 120001 : iGehrels et alJ 
120061 iLiang et al.l l2006h. but is typically neg li gible for 
short -duration GRBs (|Norris fc Bonneifeooel; lYi et all 
2006). Usually, the lag may be visualized as the time 
differences of the peaks of the "same" pulse in different 
energy bands. Statistically it can be derived through a 
cross-correlati on analysis of the p ulse profiles in different 
energy bands (jNorris et al.ll200"fj ). Technically, what is 
usually measured is the lag time At between two BATSE 
(or Swift BAT) bands E and E + AE. Mathematically, 

this corresponds to \ J^ +AE (dt/dE)dE\. It is therefore 
important to study dt/dE (or dE/dt) in theoretical mod- 
els. 

Theoretically, the leading model of the spectral lag is 
the "kinetic" effect, i.e. the delay is due to the fact that 



the observer is looking at the increasing lat itudes with 
respe ct to the line of si ght w it h time (e.g. ISalmonsonl 
2000: llokafc Nakamurall200l]: iNorris fc Bonnelll 120061: 
Shen et al.l l2005t iLu et alJ |2006[) . One can derive the 
spectral lag within this model as follows. Since the cool- 
ing time scale in the GRB emission region is very short, 
one may assume that the decay of GRB pulses is dom- 
inated by this high-latitude "curvature" effect. The co- 
moving emissivity is assumed to be uniform everywhere 
across the conical jet. Softer emission comes from higher 
latitudes (due to their smaller Doppler factor) and there- 
fore is delayed by a time t ~ (1 + z)(-Rgrb/c)(1 — cos 9) 
with respect to the emission from the line of sight, where 
9 is the angle from the line of sight. The observed pho- 
ton energy is related to the comoving one via E — T>E' ', 
where V is the Doppler factor [Eq. (J3J)] . One therefore 
has 



cE' 



dV 



dE _ 

dt Rgkb(1 + z) d(-o 
cE'P 



(1 + z) R GRB T (1- (3 cos 9f 



(11) 



The negative sign denotes "lag" , i.e. increasing E 
(harder) corresponds to a decreasing arrival time t (ear- 
lier). When 6* — > (close to the line of sight), this is 
simplified as 



dE 
~dt 



4cET 3 
(1 + z)Rgrb 
2E' 



(l + z)n 



2 E 

1 + Z T 



where 



R 



GRB 



2r 2 c 



(12) 



(13) 



is the angular spreading time, which could be related 
to the observed half width (in the decaying wing) of the 
GRB pulse 15 . Notice that -Rgrb is a value one cannot di- 
rectly measure, therefore in the above expression it needs 
to be combined with T 2 to derive r, leaving only one 
power in the T-dependence in Eq. (|12p . Another com- 
ment for this expression is that dE/dt is independent, 
i.e. for the same GRB pulse, the lag between E\ and 
(E\ + AE) should be different from that between E% 
and (E 2 + AE). This can be understood with Eq. (fTTj) 
by noticing that different E corresponds to different T>, 
and hence, different 9. For 9 > 0, dE/dt takes a differ- 
ent form than Eq. (fT2| , which is valid for the hardest (on 
axis) pulse. 

One can immediately draw the following inference from 
Eq. (fT2")) . Since dE/dt ~ E/t, one can get the lag At ~ r 
if one takes AE ~ E. This is to say, the lag time is 
comparable to the pulse width itself. So spectral lags do 
not carry direct information of the progenitor. On the 
other hand, since pulse width can be related to variabil- 
ity time scale, which may be related to the physical types 
(§3.5), one may speculate the expected spectral lags of 
the two types by the following way. Type I GRBs have 
naked central engines, so their pulse widths r are typi- 
cally much smaller than those of Type II GRBs, whose 

15 In principle, r is the upper limit of the observed half width 
in the decaying wing, since part of the tail may be buried beneath 
the next rising pulse. 
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variability time scales are longer due to the additional 
modulation of the stellar envelope. One therefore may 
expect that Type I GRBs have shorter lags than Type 
II GRBs. This is consistent with the fact that short du- 
ration GRBs (preferentially Type I) have negligible lags, 
while long duration GRBs (preferentially Type II) have 
long lags. 

Another commonly discussed argument is that Type 
I GRBs may have larger T's than Type II GRBs (due 
to a "cleaner" environment with less baryon loading), 
and that this might b e the origin of short lags (e.g. 
iNorris fc Bonnelll 12006). According to Eq. (fT2"|) , one can 



argue 



16 



dt 
dE 



oc 



1 

r 



(14) 



if different bur sts have similar r and E' (see also 
IShen et aT1l2005f ). However, the two assumptions (same 
r and E 1 ) lack physical justifications. In particular, E 1 
depends on the dissipation mechanism and the proper- 
ties (e.g. B field strength) in the dissipation region, 
which depends on the burst parameters such as L, dt, 
etc. Although the central engine variability time scales 
(r) may be arguably similar within the Type I or Type 
II category, respectively, they are considerably different 
between the two types. We therefore conclude that the 
Lorentz factor argument is not robust. Type I GRBs can 
have larger T's, but it is not the main reason for their 
short spectral lags. 

A major issue of such a kinetic (high-latitude cur- 
vature effect) model is that t he peak flux of the 
pulse drops rapidly with angle (|Fenimore et alj 119961 : 
iKumar fc Panaitescull2000T ). so that the flux is expected 
to be too low in softer bands to interpret t he observed 
flux None of the previous kinetic modelers (ISalmonson 
2000tlloka fc Nakunurall2001tlShen et alj|2005t iLu et al 



20061 ) have seriously confronted the flux predictions with 
the data (although the timing data have been well inter- 
preted by the models). One way to increase the high- 
latitude flux is to invoke a non-power-law instantaneous 
spectrum at the end of internal emission (e.g. at the 
shock crossing time in the internal shock model), e.g . 
a power law with exp onential cu toff (jZhang et al.ll2009l) 
or a Band function (|Qinl [2008) . The curvature effect 
of these models predicts that the spectral peak sweeps 
across different energy bands, making the flux not drop 
as rapidly as in the case of a power law spectrum. In- 
deed, GRB 060218 can be modele d by an evolving cutoff 
power law (|Campana et al.ll2006[ ). and GRB 050814 can 
be modeled by the curvat ure effect model inv oking a cut- 
off power law spectrum (jZhang et al.l l2009h . However, 
the light curve for a given band is always a decay func- 
tion unless the spectral index before E p is much flatter 
than -1. This is not supported by the spectral data of 
most GRBs. One is then obliged to abandon the hypoth- 
esis of a uniform jet. A structured jet with a less energy 
and/or a lower Lorentz factor at large angles from the jet 
axis and with the line of sight piercing into the wing of 
the structured jet (Zhang & Mcszaros 2002b; Rossi et al. 



16 Notice that this is different from Norris & Bonnell (2006) who 
argued At oc r — 1//2 based on the expression of the angular spread- 
ing time rather than based on the differential property dE/dt as 
discussed in this paper. 



[200l [Zhang et al.ll200l can be invoked to account for 
the observed spectral lag data. 

In reality, there might be additional mechanisms that 
are related to the observed spectral lags of GRBs, but 
the kinetic effect must exist, and may play the dominant 
role to define spectral lags in most GRBs. 

4.7. Statistical Correlations 

Observationally, some empirical correlations among 
several observed quantities have been claimed (see e.g. 
Zhang 2007 for a summary). Most of these correlations 
were discovered for long/soft GRBs. Here we discuss two 
of them that are potentially related to Type I/II diver- 
sity. 

Amati (Yonetoku) relation. Statistically, more en- 
ergetic long GRBs are harder. This is usually expressed 

in terms of E p (l + z) oc El]J so ([Amati et al.ll2002f ) and 

E p (l + z) oc L^? ao ([Yonetoku et al.l l2004h . The dis- 
persions of the correlations are large, and outliers do 
exist (e.g. the nearby GRB 980425/SN 1998bw is an 
outlier of the Amati- relation). It has been argued that 
the observe d correlations are solely due to some selec- 
tion effects (iNakar fc PirarJ l2005t iBand fc Preecel l2005t 
iButler et al J 120071) . However, the fact that the correla- 
tions are valid for most z-kn own GRBs which cove r five 
decades in isotropic energy ( Sakamo to et al.ll2006j ) sug- 
gest that there is likely underlying physics that drives 
such correlations. 

Inspecting the expressions of E p in various GRB 
prompt emission models discussed in §3.5 (see also 
Zhang & Mcszaros 2002a), one can see that although E p 
is indeed generally a function of L, it is usually also 
a function of other parameters, in particular, the bulk 
Lorentz factor T which is usually not directly measured. 
In order to interpret the Amati/ Yonetoku relations, one 
needs to introduce a rough dependence between T and L. 
We now again discuss the three prompt emission models 
and investigate how the correlations may be interpreted 
within each model. 

• For the internal shock models, the Am- 
ati/Yonetoku relations require that -f^r^ 1 is 
roughly constant. If 7 e is roughly constant, then 
the internal shock radius should be similar for 
different bursts. This suggests that T is essentially 
independent of L (since the variability time scale 
may be similar among Type II GRBs), a require- 
ment not immediately evident based on physical 
arguments. Alternatively, if "f e oc T, the relation 
can be naturally satisfied for the synchrotron 
model (k = 2). This suggests that high-r GRBs 
are more variable. In other words, while the 
maximum Lorentz factor of the outflows Tm can 
vary from burst to burst, the minimum Lorentz 
factors r m for different bursts are similar to each 
other, so that 7 e oc (Fjvf/r m ) oc F. 

• Within the photosphere model, the most relevant 
regime for Type II GRBs is the one with a stel- 
lar envelope (Eq.[5]). An interpreta tion of the 
Ama t i /Yonetoku relation can be made ([Thompson! 
120061 : [Thompson et all l2007t ) by introducing an- 
other assumption that the total energy of dif- 
ferent GRB jets is quasi-universal, a conclusion 
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reached in t he pr e-Swift era (|Frail et al.l |200H : 
iBloom et al.l 120031 ). Recent Swift observations 
suggest that the "achromatic" behavior, the de- 
manded cha racteristic of a jet break, is not com- 
monly seen (jLiang eTaT][2(¥)l . This raises the is- 
sue of interpreting some afterglow temporal breaks 
as jet breaks. Furthermore, the inferred total 
energies (after beaming correction) are found to 
a h ave larger s c atter than the pre-Swift sam- 
ple (iLiang et all l2008t iKocevski fe Butler! 120081 : 
iRacusin et al. I 120091) . In any case, if one be- 

const, using the argu- 



the pulses become na rrower with energy w i th w oc E a 
with a ~ (0.3 - 4) (iFenimore et al 1 119951: [Norris et all 



lieves E~ 



E 9 2 



ment that baryon sheath near the breakout ra- 
dius leads to T„ ~ fl" 1 oc E^J so (|Thompsonl 120061 
iThompson et alj [2007) . one can translate Eq.© 
into E p oc LljJ so oc E^J so by taking the triv- 
ial proportionality L oc L 7) j s0 oc E. 



The 

first proportionality is based on the fact that 
the GRB efficiency is not a function of L 7 ,j s0 
(jLloyd-Ronning fc Zha ng 2004); while the second 
proportionality is based on the fact that the dis- 
persion of T90 is not large and that Tqq is not cor- 
related to L l iso . In view of the fact that GRB 
080916C disfavors th e photosphere origin of E p 
(Zha ng fc Pe'erl I2009D , we regard this interpreta- 
tion as no longer attractive. 

• The magnetic dissipation model lacks a robust pre- 
diction for E p . In any case, similar to the other 
two models, the Amati/Yonetoku relation can be 
satisfied if one assigns a particular T — L correla- 
tion. For example, the specific model described in 
Eq.([T0"]) requires T oc L 1 / 2 . 

When the optical afterglow light curve temporal break 
(topt) is included, a tighter correlation involving E v and 
Ei.is o is obtained (|Ghirlanda et al. 2004: Liang & Zhang 
200o| for some long duration GRBs. However, the phys- 
ical connections between the prompt emission properties 
(E p and E lyiso ) and the afterglow properties (t opt ) or the 
global collimation degree of the jet are not straightfor- 
wardly expected. Furthermore, optical afterglow tempo- 
ral break data of Type I GRBs are still poor to draw any 
conclusion. We therefore do not discuss these relations 
in this paper. 

Luminosity-lag relation. INorris et alj |2000) dis- 
covered a relation between the gamma-ray peak luminos- 
ity L^ iso and the spectral lag At between the BATSE 
Channel 1 (25-50 keV) and Channel 3 (100-300 keV), 
i.e. iso oc (At) -1 . The relation was refined by 
iGehrels et al.1 (|2006l ) who corrected the observed spec- 
tral lag to that between two common bands in the cos- 
mic proper rest frame of GRBs. This correction takes 
two steps. First, the observed spectral lag At between 
the bands from E to (E + AE) can be expressed as 
At = At rest (E rest )(l + z), where At rest (E rest ) is the 
spectral lag between E rest to (E rest + AE rest ) in the 
cosmic proper rest frame. Second, what one cares about 
is the intrinsic spectral lag between a common rest frame 
energy interval, e.g. from E to (E + AE). One needs an 
additional relation between At res t(E) and At rest (E rest ). 
There is no universal relation for this, but there is an 
empirical relation between pulse width and energy, i.e. 



120051: ILiang et alj 120061: IZhang fc Qinll2008l ). Assuming 
that the spectral lag is pr oportional to the pulse width 
(see e.g. lNorris et al.ll2005h . one may derive At rest (E) = 
At res t(E rest )(E j E rest )~ a = At res t(E rest )(l + z) a . This 
finally gives 

At rest (E) = At(E)(l + z) a-1 . (15) 

In IGehrels etaLl ((2006ft . a ~ 1/3 was adopted, and a 
correlation 



L? lso cx{At(l + zy 



At 



(1 + z) 2 / 3 



(16) 



is found for a sample of long GRBs, with S ~ 1. Outliers 
do exist for long GRBs, and short GRBs are noticeably 
off the track of the correlation. 

Is there an underlying physical mechanism that jus- 
tifies the observed iso oc (At rest )~ s correlation? It 
is not obvious based on the theoretical arguments above. 
According to Eq. (fT2"]) , the intrinsic lag is related to the in- 
trinsic variability time scale of the burst. So a L— lag neg- 
ative relation may be related to another, probably more 
intrinsic L — V positive relation, where V is the variabil- 
ity parameter. Technically, th ere are different definitions 
of the variability parameter (IFenimore fc Ramirez-Rur3 
I2000t iReichart et all 120011: iGuidorzi et all 120061 ). but in 
any case, a more variable light curve (high V) would 
have shorter variability time scales (corresponding to t), 
and hence, shorter spectral lags (At). Observationally, 
indeed a positive L — V relation is observed, although 
with a large scatter (IFenimore fe Ramirez-Ruizj 120001 : 
IReichart et aLll200H IGuidorzi et al.ll2006l) . The interpre- 
tation of this correlation within the internal shock model 
invokes several assumptions: (1) The smallest variabil- 
ity time scale is defined by the collisions above the pair 
photosphere: (2) The true jet energy is quasi-universal 



( Frail et all I200l1 ) ; (3) Narrower jets have h igher T's 
( Kobavashi et al.l l2002t iMeszaros etTafl [2002 i 17 . This 
is a relevant interpretation to the observed iso oc 

(At r est)~ S relation. However, in view of the assump- 
tions invoked in the reasoning (the above three as well as 
the assumption that the w — E correlation is similar to 
At — E correlation as invoked earlier), we expect that the 
correlation should not be very tight, and may not follow 
the same simple power law. This is consistent with the 
data (see §5 for details). 

Would Type I GRBs satisfy a similar correlation? In 
principle one can expect so if the set of assumptions 
discussed above are satisfied. In reality, Type I GRBs 
may not effectively develop a pair photosphere, both be- 
cause of their preferred higher Lorentz factors (clean en- 
vironment) and because of their lower total energy. This 
would make the observed variability time scale trace the 
central engine variability time scale, which would not 
vary significantly among bursts. 

Another model to interpret the iso oc (At rest )~ s 
correlation in vokes varying Doppler factors among dif - 
ferent bursts (|Salmonsonll2000l : lloka fc NakamurallioOll) . 

17 The third assumption may be in conflict with the explanation 
of the Amati/Yonetoku relation within the standard internal shock 
model as discussed in §4.5. To interpret that relation, one requires 
no dependence of T on L (and hence, on 8). 
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These models assume a universal comoving properties 
of ALL GRBs, and invoke the off-beam geometry to in- 
terpret longer durations and spectral lags. As already 
discussed, the comoving properties of GRBs depend on 
many parameters. The off-beam model is not supported 
by early afterglow observations. We therefore regard 
these early models invoking pure geometrical effects as 
no longer favorable in view of the recent observational 
progress. 

4.8. Energetics and Beaming 

Type II GRBs are generally expected to be more en- 
ergetic than Type I GRBs. In the standard BH-torus 
central engine model, the total energy of the burst is 
positively correlated with the total available fuel in the 
torus. Massive stars are much more abundant in mass, 
which can reach ~ 10M Q of fuel in total. On the other 
hand, a NS-NS merger system has a total energy budget 
of ~ 2.8M©. After the prompt collapse, the available 
fuel in the torus is of order ~ 0.1M©. A BH-NS merger 
system has even less fuel to begin with (~ 1.4M©). It 
would reach a similar total energy budget in accretion as 
the NS-NS system. In these models, Type I GRBs are 
expected to be 10-100 times less energetic than Type II 
GRBs. Alternatively, GRBs may be powered by the spin 
energy of the central object (a rapidly spinning BH or 
NS). For the case of a NS central engine, Type I GRBs 
may reach similar energies as Type II GRBs. For the 
case of a BH engi ne whose spin energy is ext rac ted vi a a 
magnetic torque (|Blandford fe Znaieklll977t [L1 I2002T ). a 
Type II GRB is again expected to be ~ 10 times more 
energetic than a Type I GRB if it is powered by a NS- 
NS merger, again because of the more massive BH in the 
Type II GRB. A BH-NS merger Type I GRB, on the 
other hand, may reach the same energetics as Type II 
GRBs if the initial BH in the binary system is massive 
enough and has a large enough angular momentum. 

In order to relate this theoretically motivated total en- 
ergy budget to the observed energy, one needs to in- 
troduce the beaming factor. The standard GRB jet 
models in yoke a conical jet with uniform energy dis- 
tribu tion (|Meszaros et al.lll998t lRhoaddll999t ISari et all 
1999). More complicated (maybe more realistic) jet mod- 
els invo ke distributions of jet energy with angle from the 
net axis (iMeszaros et aDll998t IZhang fc Meszarosll2002bT : 
iRossi et al.l l2002t IZhang et al l 120041) . Theoretically, it 
is difficult to model how jets are launched from the 
central engine. On the other hand, one can speculate 
about the collimation angle of jets from two types of 
GRBs from the theoretical point of view: Type II GRBs 
should tend to have narrower jets than Type I GRBs 
due t o the additional collimation of the stellar enve- 
lope (IZhang et al.ll2003D Type I GRB jets tend to be 
broader ( Alov et al.ll2005l ). Observationally this predic- 
tion has not been tested statistically. Observations of 
some individual bursts seem to support this picture. For 
example, the Type I Gold Sample burst GRB 050724 
was found to have a beaming angle wi der than ~ 25° 
(|Grupe et alJl2006b iMalesani et al.ll2007h . Type II GRBs 
on the other hand, have a typical beaming angle of ~ 5° 
(iFrail et all 120011: [Bloom et all l2^0ljLiang_ etlll 120081 : 
iKocevski fe Butlerl 120081: iRacusin et al]|2009t ). Opposite 
cases are also observed in some GRBs. For example, the 
short GRB 051221 (a Type I candidate) has a narrow 



jet with 9j ~ 4° - 8° (|Burrows et al.ll2006D . This is con- 
trary to the theoretical expectation if it is indeed a Type 
I GRB. On the other ha nd, the Type II GRB 060729 
(|Gruoe et al.ll2007L l2009ah may have a large opening an- 
gle since its X-ray afterglow keeps decaying without a 
break for hundreds of days. 

The observed "isotropic" energy is the total energy di- 
vided by the beaming factor (2 • -kO^/A-k = 6^/2 for uni- 
form jets under the small angle approximation). Con- 
sidering a factor of ~ (20/5) 2 ~ 15 difference in the 
beaming factor, the isotropic energy of a typical Type I 
GRB should be a factor ~ 100 — 1000 times lower than 
that of a typical Type II GRB. If a Type I GRB has a 
high isotropic luminosity/energy comparable to that of a 
typical Type II GRB, one must then demand a very nar- 
row jet, with an opening angle even smaller than that 
of Type II GRBs. A BH-NS central engine with the 
Blandford-Znajck mechanism as the engine power can 
ease the constraint, but a comparable beaming angle to 
Type II GRBs is nonetheless needed. 

4.9. Afterglow Properties 

Broadband afterglow emission has been interpreted as 
the external forward shock emission as the fireball is de- 
celerated by the circumburst medium. An ideal observa- 
tional campaign can lead t o diagnostics of the c i rcum- 
burst medium properties (|Panaitescu fe Kumarl l2002t 

I Yost et al.|[2003h . which c an shed light on the progenitor 
system of the GRB (e g. iFan et alJ l2005t iGreiner et all 
I2009al IXu et all 120091 ). In particular if a stratified 
stellar-wind-type me dium (n oc R~ 2 ) (|Dai fc Lul 119981 : 
IChevalier fc Lill20"00[ ) is identified, the burst can be iden- 
tified as a Type II GRB. The case of a constant density 
medium is, however, less informative. Although Type I 
GRBs are expected to reside in such a medium, some 
Gold-Sample Type II GRBs have been found to reside in 
a constant medi um as well (jPanaitescu fc Kumarl 120021 : 
lYost et all 120031) . The mechanism of forming such a 
medium before the death of the massive star is unknown. 
In any case, one needs more information to establish the 
association of a burst with a physical model category 
if it goes off in a constant density medium. For exam- 
ple, the afterglow luminosity of Type I GRBs should be 
systematically lower than that of Type II GRBs due to 
the expectation of both a lower medium density in the 
merger environment (relevant for v < v c ) and a system- 
atically lower blastwave energy (jPanaitescu et al.l 120011 : 
IFan et alJl2005HKann et"aT1l2008h . A Type I GRB can be 
even "naked" (i.e. no detectable afterglow) if the ambient 
medium density is low enou gh. Such GRBs are i ndeed 
observed (e.g. GRB 051210. lLaParola et al.|[2006l ). 

A related issue is the GRB radiative efficiency. Ob- 
servations and theoretical modeling suggest that the ef- 
ficiency is similar for both Type I and Type II GRB s 
(|Zhang et all l2007at iBergerl 120071 iGehrels et all 120081) . 
so that it cannot be regarded as a useful criterion to tell 
the model category that is associated with a particular 
GRB. 

4.10. Redshift Distribution 

Statistically, redshift distributions of Type I and Type 

II GRBs should be different. Type II GRBs generally 
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trace the star- forming history of the universe . Type I 
GRBs are expected to be "delayed" with respect to star 
formation due to the long merger time scale asociated 
with the shrinking of the binary orbits due to gravita- 
tional radiation (jBelczynski et al.ll2006l ). On average, it 
is expected that the mean redshift of Type I GRBs is 
lower than that of Type II GRBs. 

4.11. Luminosity Function 

The luminosity function of long duration GRBs is 
categorized by a broken power law with a break 
~ 10 52 erg s" 1 dGuetta et al.l l2005t iLiang et all 120071 : 
I Virgili et al.ll2009afh Below the break the power law in- 
dex is > — 1 , while above the break the power law index is 
< — 2. Low luminosity GR Bs may form a d istinct bump 
at L < 1 49 er g s _1 (|Liang et all 120071 : IVirgili et all 
I2009al lDa1l2009D . As argued below, most long GRBs 
are Type II GRBs, so this luminosity function may be 
regarded as that of Type II GRBs. There is however no 
direct theoretical reason for such a luminosity function. 
For Type I GRBs, the luminosity function has not been 
studied in detail due to the li mited sample with r edshift 
measurements so far (but see IVirgili et al.1 120 09b'). The 
luminosity funct ion of the B ATSE shor t /hard GRB sam- 
ple was studied (jGuetta fc PiranlfeOOGf ). but as discussed 
below, it is not justified that this population is identi- 
cal to the Type I population. This issue will be further 
discussed in §5.4 and §6 below. 

4.12. Gravitation Wave Signals 

Probably the most definite criterion to differentiate 
Type I GRBs from Type II GRBs is through detect- 
ing their gravitational wave (GW) signals. Although 
the GW signature of a Type II GRB is highly uncer- 
tain (e.g. iKobavashi fc Meszarosll2003h . the wave forms 
of NS-NS and N S-BH mergers are well predicted (e.g. 
iDalal et al.ll2006l ). Detections of these signals would un- 
ambiguously associate some GRBs to the Type I model 
category. However, this criterion can only be applied 
in the future when the GW detectors reach the desired 
sensitivities. 

5. TYPE I AND TYPE II SAMPLES AND THEIR 
STATISTICAL PROPERTIES 

5.1. Sample selection 

The above consideration suggests that theoretically 
there is no handy, distinct criterion that can be used to 
immediately determine the physical model category that 
a burst is associated with. In this section, we attempt 
to explore the topic further from the observational point 
of view. The current standard approach is to use three 
criteria, "duration" , "hardness" , and (when available) 
"spectral lag", to categorize bursts as "long" (implic- 
itly assumed to be associated with "Type II" ) or "short" 
(implicitly assumed to be associated with "Type I" ) , and 
use these samples to explore the sta tistical prop e rties of 
other observational properties (e.g. iNakarl l2007t iBergerl 
2009). In other words, it is often implicitly assumed that 

Long/soft/long lag = Type II 

Short/hard/short lag = Type I. (17) 

18 Metallicity may play additional role to select Ty pe II G RBs 
(e.g. I Wolf fc Podsiadlowskil [20071: INuzaet alj|2007t H^f l2008h . but 
the issue is inconclusive. 



The problem with such an approach is that these crite- 
ria may not be always reliable. A notable example was 
GRB 060614, which is a long GRB but is very likely as- 
sociate d with Type I (jGehrels et al.ll2006t I Zhang et all 
l2007bh . Another complication is related to GRB 080913, 
GRB 090423, and som e other intrinsically short GRBs 
(e.g. iLevan et al.l I2007D . These GRBs can be detected 
as short/hard GRBs if their rcdshifts were low enough, 
but their physical properties are more close to those of 
Type II GRBs. Some short/hard GRBs (e.g. 060121, de 
Ugarte Postigo et al. 2006) are very energetic, which are 
not easy to be accommodated within the Type I progen- 
itor models. 

In this paper, we adopt an alternative approach. In- 
stead of sticking to the observed gamma-ray properties, 
we adopt the observational criteria that are directly re- 
lated to the progenitor systems to select the samples, and 
then go back to investigate the other properties (includ- 
ing duration, hardness, spectral lag, etc) of the samples. 
The advantage of this approach is that we can start with 
those GRBs whose progenitor systems are more confi- 
dently inferred. We can then use them to verify whether 
the ansatz Eq. (p^Tjl is justified. 

We define the following three samples based on the 
criteria detailed below. 

Type II Gold Sample. This sample is defined such 
that at least one of the following two criteria are satisi- 
fied. 

1 . There is a spectrally confirmed SN association with 
the GRB; 

2. The specific star forming rate (SSFR) is very high 
(to be specific, the SSFR satisfies log SSFR > -0.2 
or SSFR> 0.63 Gyr -1 in the sample of Savaglio et 
al. 2009); the GRB location does not have a large 
offset from the center; and there is no stringent 
upper limit on the existence of a SN associated with 
the GRB. 

Notice that the GRB properties (duration, hardness and 
lag) are not the considerations to define the sample. 
Since not many GRBs have host SSFR information pub- 
lished, this sample is by no means complete, and there 
should be many more Type II GRBs that are not in- 
cluded. The purpose of selecting this sample is to use 
the most stringent criteria to investigate how the best 
Type II GRB candidates look like. As a result, we do 
not include the GRBs that have a claimed SN bump 
in the optical light curve but no confirmed SN spec- 
troscopic signature. The threshold of SSFR is arbi- 
trary. This lim i ting v alue was chosen because Table 1 1 of 
iSavaglio et al.l (|2009f l has a mix of long and short GRBs 
for log SSFR(Gpc _1 ) < —0.3, which is the regime where 
confusion arises. The lower bound log SSFR > —0.2 
can be regarded as a safe line above which GRB hosts 
have very active star formation. One exception is the 
short duration GRB 051221A (jSoderberg et all 120061 : 
iBurrows et all l200l . The SSFR value (log SSFR > 
0.804) is way above the threshold. However, since deep 
sear ches have ruled out th e association of a 1998bw-likc 
SN dSoderberg et al.ll200fih . we do not include it in the 
Type II sample, and will include it in the "Other short 
hard sample". We note that many Swift long GRBs 
should be associated with Type II. However, since no 



Physical origins of GRBs 



15 



published SSFRs are available for most of them, we re- 
frain from including them in the Type II Gold Sample 19 . 
This sample should be expanded significantly later when 
the host galaxy information of the Swift GRBs is re- 
leased. Right now the Type II Gold Sample includes 
33 GRBs (Table 1 Top Panel). This is already a large 
enough sample to study the statistical properties of Type 
II GRBs. 

Type I Gold Sample. The Gold Sample of Type 
I GRBs is defined by at least one of the following two 
criteria. 

1. The host galaxy is elliptical or early type; 

2. The GRB location has a relatively low local SSFR, 
or a large offset from the center of the host galaxy; 
and deep searches reveal stringent upper limits on 
the existence of an underlying SN. 

Again the GRB properties (dura tion, hardness, lag) are 
not considered. Some argum ents (jBelczynski et al.ll2006t 
IZheng fc Ramirez-Ruij[2Q07l ) have suggested that a frac- 
tion of Type I GRBs may be located in star forming re- 
gions of star forming galaxies. Our criteria do not select 
those, since we do not demand completeness of sample se- 
lection. After systematically checking the archival data, 
we only identify 5 bursts in the Type I Gold Sample: 
GRBs 050509B, 050709, 050724, 060614 20 , and 061006 
(Table 1 Middle Panel). The details of individual GRBs 
are presented in the Appendix. 

Other SGRB Sample. Most short/hard GRBs in 
the Swift era satisfy neither of the two criteria of the 
Type I Gold Sample. Some of them do not have their 
host galaxies convincingly identified. Others have host 
galaxies with active star formation. These GRBs are 
usually regarded as Type I candidates simply because 
they are "short/hard" . There could be a good fraction 
of Type I GRBs in this sample, but we are not sure that 
they can ALL be associated with Type I. Since we define 
the Gold Samples not based on the GRB properties, we 
leave these bursts in a separate sample, without speci- 
fying whether they are associated with Type I or Type 
II. There are 20 bursts in this sample (Table 1 Bottom 
Panel) . The details of individual GRBs are presented in 
the Appendix. 

5.2. Duration-Hardness Distribution 

Figure presents the traditional T 90 -hardness ratio 
(HR) plot of GRBs. Superimposed on the BATSE data 
(orange small dots) are the three samples defined above: 
Type II Gold Sample (blue), Type I Gold Sample (red), 
and other SGRB sample (green). The HR is defined as 
the fluence ratio between (50-100) keV and (25-50) keV. 
For BATSE bursts, this corresponds to the fluence ratio 

19 Besides those included in the ISavaglio et al. (2009) sample 
(which covers from GRB 970228 to GRB 061126), we only in- 
clude GRB 080520 and GRB 060602A based on the SFR criterion. 
They have a high SFR (though SSFR is not measured) typical 
for other Type II Gold sample G RB host galaxies. Fo r example, 
GRB 080520 has ~ 15M Q vr" 1 UMalesani et al l 120081) , which is 
comparable to the highest in the Savaglio et al. (2009) sample. 

20 In the literature GRB 060614 is usually taken as a contro- 
versial candidate for Type I. This was mainly because of its long 
duration. We do not consider duration as a criterion when select- 
ing the Gold Sample. This burst satisfies the criterion #2 of the 
Type I Gold Sample. 



between channel 2 and channel 1. For other detectors 
(HETE-2, Swift/BAT, Konus/Wind, INTEGRAL) with 
different detector energy bands, we perform spectral fits 
and use the fitted model to derive the HR. Besides the 
observed points (open symbols), we also plot the cor- 
responding "rest-frame" points (filled symbols) for each 
burst. The HR is then defined as the flux ratio between 
the rest-frame (50-100) keV band and (25-50) keV bands, 
which is again derived from spectral fitting. For a power 
law fit, the rest frame HR is the same as the observed 
one. For a curved spectrum (e.g. a Band function or an 
exponential cutoff power law), the two can be different. 
The Tg values are energy- and detector-dependent. We 
do not make efforts to convert all Tgo to the BATSE- 
band, since this requires time-dependent spectral anal- 
yses and extrapolations, and for many bursts the data 
quality is not sufficient to perform such an analysis. In- 
stead we simply plot Tgo measured by different detectors 
(e.g. Swift and HETE). The correction to the BATSE- 
band Tgo is usually not significant for most long GRBs, 
but could be significant to those GRBs with soft ex- 
tended emission. Traditionally, the "rest frame" Tgo are 
not used to defined long vs. short for a particular GRB. 
We present them here just to show how the intrinsic dis- 
tribution may differ from the observed one. To derive the 
rest- frame Tgo st , we simply divide the observed value by 
(1 + z). More rigorously one needs to again take into ac- 
count the light curve evolution with energy. This again 
requires a time-dependent spectral analysis. Since most 
bursts do not have such detailed information, and since 
the correction would not be significant for most bursts, 
we neglect this correction for the sake of simplicity and 
uniformity. For short GRBs with extended emission, we 
use circles to denote the short spikes only (excluding the 
extended emission) , while using squares to denote the full 
emission with extended emission included. These two lo- 
cations for the same burst with and without extended 
emission are connected by lines. Since the mean HR is 
derived, the HRs including extended emission are usually 
smaller than those without, as the extended emission is 
typically softer than the initial short spikes. 

From Fig[2]one can make the following interesting ob- 
servations. First, the Type II GRBs are generally long, 
and they well represent the long/soft population of the 
BATSE GRBs in the T 90 -HR plane. However, some Type 
II GRBs have a duration close to the 2-second separation 
line, and their intrinsic duration can be shorter than 2 s 
(e.g. GRB 040924 with T 90 = 2.39 ±0.24 s at z = 0.858, 
and GRB 080520 with T 90 = 2.82 ± 0.67 at z = 1.545). 
iLevan et al.l (|2007l ) also discussed a sample of apparently- 
long, intrinsically-short GRBs. Secondly, four out of five 
Type I Gold Sample GRBs are not strictly "short". Ex- 
cept GRB 050509B, all the others have extended emission 
aside from the initial "short/hard" spike. The spike it- 
self is longer than 2 s for GRB 050724 and GRB 060614. 
All 5 Type-I Gold Sample bursts have a moderate HR. 
None has an extremely hard spectrum. Thirdly, the 
Other SGRB Sample fills in the short/hard region in the 
Tgo— HR diagram more uniformly, suggesting that it rep- 
resents the BATSE short/hard sample well. Some bursts 
in the sample also have extended emission. 

5.3. Empirical correlations 
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Figured displays the E p — E^s so (Amati) relation of 
the three samples. The spectral parameters are collected 
from the published papers or GCN circular reports (see 
Table 1 for references). For those GRBs with extended 
emission (including Type I Gold Sample GRBs 050724, 
060614, and 061006), we only consider the short hard 
spikes. For all the bursts, the isotropic gamma-ray en- 
ergy (E 7y i SO ) is calculated in the GRB rest-frame 1 — 10 4 
keV band through extrapolation based on the spectral 
parameters. We can see that most GRBs in the Type II 

1 /2 

Gold Sample indeed follow the E p cx E^' iso (Amati) rela- 
tion. However, there are three noticeable outliers: GRB 
980425, GRB 031203, and GRB 050826. The first two 
are nearby low- luminosity (LL) GRBs, whic h have been 
argued to be from a distin c t populati on (e.g. lLiang et al.1 
120071 IVirgili et al.l l2009at [Dai 2009). Another nearb; 



LL GRB 060218 is a soft burst ()Campana et al.ll200 
and satisfies the Amati-relation well. GRB 050826 with 
Tqo ~ 35 s is an intermediate Type II GRB between 
the more "classica l" Type II and the nearby LL-GRBs 
(|Kann et al.l l2007h . and deviates from the relation. We 
also pay special attention to the two intrinsically short 
Type II GRBs. Although GRB 040924 is right on the 
Amati-relation track, GRB 080520 seems to be slightly 
off the track. The Type I Gold Sample and the Other 
SGRB Sample are populated above the conventional 
Amati-relation track. Since many short/hard GRBs have 
E p outside the BAT band, their E p error b ars are large. 
The y alues in our analyses are adopted from lButler et all 
(2007). In any case, it seems that they follow a sepa- 
rate track with a shallower slope than the Amati-relation. 
Excluding GRBs 080913, 090423 and 060121 (which are 
likely Type II, see §6.2), a best fit to the Type I Gold 
and Other SGRB samples lead to a slope 0.34, with the 
3(7 limits of the slope as (0.15-0.53) (see FiglSi). GRB 
080913 is marginally within the 3<r regions for the Type II 
Amati-relation, but is also consistent with this new track 
defined by Type I and other short/hard GRBs within 3a. 
GRB 090423 ali gns with the Typ e II Amati-relation more 
closely (see also lLin et al.1120091 ). 

A likely reason that the Type I and the Other SGRB 
Samples deviate from the Amati relation of Type II 
GRBs is simply because they have shorter durations 
so that they have smaller I? 7i i SO values than the Type 
II GRBs with a similar E p . To test this, we plot the 
E p — iso relation (Yonetoku relation) in Fig[3b. We 
can see that the distinction between Type II and Type 
I GRBs becomes less significant, although the corre- 
lation now has a much larger scatter. Noticing the 
large error bars of the Type I and Other SGRB Sam- 
ples, one may conclude that there is no distinct differ- 
ence among the three samples as far as the Yonetoku 
rel ation is concerned. A similar conclusion was drawn 
by iGhirlanda et all (|2009l ) in an anaylsis of the BATSE 
GRBs. 

Figure [4^i displays the luminosity-spectral lag dia- 
gram of GRBs with the three samples plotted. A 
group of Gold Sample Type II GR Bs indeed define a 
lA.iso oc (At re st )~ s correlation track (|Norris et al.ll2000l : 
iGehrels et al.1 |2006[ L although several low- luminosity, 
long-lag GRBs lie below the extrapolation of th e track 
(see also IGehrels et all 120061: iLiang et al.l l2006t) . Gold 
Sample Type I GRBs are clustered at the lower left cor- 
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Fig. 3. — (a) The E p — E lt i ao diagram of the three samples of 
GRB discussed in the paper: Type II Gold Sample (blue), Type I 
Gold Sample (red), and other short/hard GRBs (green). Two pos- 
sible redshifts z = 4.6, 1.7 for the short GRB 060121 are adopted, 
which satisfies the relation well (unlike other short/hard GRBs). 
GRB 080913 and GRB 090423 (cyan) are also plotted for compar- 
ison. The best-fit E v — E lt i ao correlations for both Type II and 
Type I/Other SGRB samples are plotted (solid lines) with the 3<r 
boundary (dashed line) marked, (b) The E v — igo diagram. The 
same convention has been used. 



ner. This is as expected: short durations define short 
lags, and smaller energy budgets define lower luminosi- 
ties. About half of the "Other SGRBs" are clustered 
close to the Type I Gold Sample, suggesting that they 
may be associated with Type I as well. Some others fill in 
the gap between the Type I and Type II Gold Samples. 
In particular, GRB 060121 lies r ight on the track for both 
putat ive redshifts 1.7 and 4.6 (|de Ugarte Postigo et al.1 
I2006f) . GRB 070714B is also close to the track. The 
SN-less GRB 060505 clusters with other nearby low- 
luminosity Type II GRBs. Finally, the two high-z GRBs 
080913 (notice that only the upper limit of spectral lag 
is derived) and 090423 are consistent with satisfying the 
, so — lag correlation of Type II, but are also consistent 
with the zero-lag trend of Type I/Other SGRB. 
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Fig. 4. — (a) The L 1 ^ iao — lag diagram. Same convention as 

Fig.3 is adopted. GRB 080913 and GRB 090423 satisfy both the 
correlation defined by Type II GRBs and the "zero lag" trend de- 
fined by Type I and Other SGRB Samples. Two possible redshifts 
z = 4.6, 1.7 for the short GRB 060121 are adopted, which satis- 
fies the correlation well (unlike other short/hard GRBs). (b) the 
lag — T90 (intrinsic) diagram of the three samples. The same GRBs 
with/without extended emission is connected by dotted lines. The 
spectral lags of these GRBs are for the short/hard spikes only. A 
positive correlation between duration and spectral lag is derived 
(dashed line). See text for details. 

As discussed in §4.7, the luminosity lag relation may 
be related to the variability-luminosity relation, and may 
be more relevant to Type II GRBs. On the other hand, 
the physical origin of the relation is not clearly under- 
stood and is based on many assumptions. Although the 
correlation may be taken as a reference, it may not be 
taken as the definite criterion for judging the physical 
origin of a GRB. 

Based on the high-latitude-effect interpretation of 
spectral lag (§4.7), one expects that short spectral lags 
should be related to short angular spreading times. The 
latter corresponds to the width of individual pulses. If 
the number of pulses do not fluctuate significantly among 
bursts, one would also expect a rough correlation be- 
tween spectral lags and durations. In Fig (4Jd we display 



the T 90 /(l + z) - lag/(l + z) 2 / 3 dia gram of the three 
samples of bursts. Again points of the same burst with 
and without extended emission are connected by lines. 
We investigate a possible correlation between duration 
and spectral lag. Since the spectral lags are defined for 
the short /hard spikes only for those GRBs with extended 
emission, we use Tqq excluding the extended emission for 
those bursts. A positive correlation between T90 and lag 
with slope 0.94 ± 0.14 is obtained, with the Spearman's 
rank correlation coefficient r = 0.735, corresponding to 
a chance probability P < 10~ 4 . This is consistent with 
our naive expectation, suggesting that spectral lags are 
closely related to durations, and may not carry additional 
information in defining the categories of GRBs. 

5.4. Luminosity and Redshift Distributions 

Figure Eta) and (b) display the observed 2-dimensional 
luminosity-redshift (L^ iso — z) and energy-redshift 
(E 1 ,i SO — z) distributions of the three samples. GRBs 
in the Type I Gold Sample are all at z < 0.5. Includ- 
ing the Other SGRB Sample, the upper boundary of z 
reaches - 1 (except GRB 060121). The Type II GRBs 
have a wider span of redshift distribution, with the peak 
around z ~ 1. In terms of luminosity distribution, the 
Type II GRBs on average are ~ 2 orders of magnitude 
more luminous than the Type I GRBs. Type I GRBs can 
at least reach a luminosity of I? iso ~ 2.5 x 10 51 erg s _1 
(for the Type I Gold GRB 061006). Including the Other 
SGRB Sample, several short GRBs (070714B, proba- 
bly 060313, and especially the latest GRB 090510) can 
reach iso ~ 10 52 erg s _1 . GRB 060121 even reaches 

^7 iso ~ 10 53 — 10 54 er S s_1 f° r ^e two fiducial redshifts 
in discussion. This luminosity is high even for Type II 
GRBs. GRB 080913 has L P l lso ~ 1.2 x 10 53 — 



o3 



erg s 
ava et al.l 



GRB 090423 has L? lso ~ 1.88x10"" erg s 
2009). Both are moderate to high luminosities for Type 
II GRBs, and are very high when compared with the 
Type I and Other SGRB Samples (except for GRB 
060121). In the E 1 ^ so — z diagram, the separation be- 
tween Type II and Type I is more distinct, with most 
SGRB sample bursts lying below the Type II distribu- 
tion. But GRB 080913 and GRB 090423 become moder- 
ate in the Type II Sample due to their intrinsically short 
durations. The clearer separation between Type II and 
Type I/Other SGRB Samples is mainly due to the short 
duration of the SGRB sample, which makes them less 
energetic. However, GRB 060121 is still as energetic as 
the average Type II GRBs. 

5.5. Afterglow Properties 

Figures [5] and [7] present the intrinsic afterglow light 
curves in the X-ray and optical bands for the three sam- 
ples. Figure [6] presents the rest-frame 2 keV specific lu- 
minosity light curves. Since many Type II Gold Sam- 
ple GRBs are pre-Swift, we do not have many Type 
II X-ray light curves. The ones that are plotted in- 
clude two low luminosity GRBs (060218 and 050826) and 
two intermediate-to-high luminosity GRBs (080520 and 
050525 A). These do not fully represent the Type II GRB 
X-ray afterglow properties. In order to compensate for 
this weakness of sample selection, we also overplot the 
X-ray light curves of a group of early Swift long GRBs 
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Fig. 5. — (a) The jgo — z diagram, and (b) the E 1 ^ so — z 
diagram of the three samples. The same convention as Fig. 3 is 
adopted. 

in the sample of iNousek et a l. (2006). Since we already 
demonstrated that the Type II Gold Sample represents 
the BATSE long GRBs well, we assume that the Nousek 
Sample represents the Type II GRB X-ray afterglows 
well. We can see that these bursts occupy the upper 
portion of the light curve space in FigJSJ By contrast, 
the Type I Gold Sample occupy the lower portion, and 
the Other SGRB Sample populate in between with much 
overlap with both Gold Samples. Low luminosity Type 
II GRBs have luminosities comparable to Type I Gold 
Sample GRBs. 

Figure [7| presents the optical light curves with cor- 
rected i? e -magnitude by moving all GRBs to z = 1 
(|Kann et alj 120071 120081 ). One big difference between 
these optical light curves and the X-ray light curves 
(Fig. [6]) is that most Type II GRBs are represented, 
exceptions being those GRBs that had negligible opti- 
cal afterglows but strong supernovae signatures (GRBs 
980425, 031203, and XRF 060218), dark GRBs, where 
the optical emission was probably totally supressed by 
line-of-sight extinction in the host galaxy (GRBs 990506, 
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t (seconds after trigger in the observer frame assuming z^l ) 

Fig. 6. — The rest frame 2 keV X-ray afterglow luminosity light 
curves of GRB 0809f3, GRB 090423, and the three samples. All 
bursts are placed at z = 1. The color scheme is the same as in 
the other figures. Since most Type II Gold Sample bursts are pre- 
Swift ones and have no X-ra y light curves, we als o add the 2-known 
long GRBs in the sample of Nousek et al. (2006) (grey), which are 
generally believed to be Type II GRBs. GRB 080913 and GRB 
090423 (cyan) both have bright X-ray afterglows typical of Type 
II GRBs. 




t (days after burst in the observer frame assuming z = 1 ) 



Fig. 7.— The rest frame optical light curves of GRB 080913, 
GRB 090423, and the three sample s. The color scheme is th e same 
as in the other figures. Similar to IKann et~aTl ((20071 . [2001) . they 
are plotted at a common redshift of z = 1. As with the X-ray 
light curves (Fig. [6jl, the optical afterglows of the Type II Gold 
Sample GRBs are clearly more luminous than those of the Type 
I Gold Sample and the Other Short-Hard Sample. The latter two 
populations are in good agreement with each other. GRB 060121 
is the single short-hard GRB which is optically highly luminous. 
GRB 080913 and GRB 090423 both have bright optical afterglows 
typical of Type II GRBs. 

000210, 020819B, 051022), and some with very sparse op- 
tical data (XRF 020903, GRBs 030528, 05082 6, 060602A 
080520). M ost data have been taken from IKann et al.l 
(|2007l 120081 ), where the methods of creating the intrin- 
sic light curves are also presented. Similar to the X-ray 
light curves, the Type II GRB afterglows form a much 
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more luminous gro up than the Type I GRB afterglows 
(|Kann et al.ll2008l ). The light curves of Type I Gold 
GRBs and those of most Other SGRBs overlap, indicat- 
ing that they are likely drawn from the same population. 
The most prominent exception is ag ain GRB 060121 w ith 
an optically luminous afterglow (see lKann et al.ll2008l for 
more details), which is comparable to the afterglows of 
Type II GRBs. 

For both X-ray and optical afterglows, GRB 080913 
and GRB 090423 have a luminosity comparable to or 
higher than the average luminos i ty of the Type II GRB 
afterglows (iGreiner et alll2009al ISalvaterra et alJ 120091 : 
iTanvir et alJl2009t ). 

6. NATURE OF SHORT /HARD GRBS 

Based on the theoretical considerations and the statis- 
tical analyses presented above, in this section we attempt 
to address the question whether the ansatz Eq. (fTT|) is 
valid, i.e. whether Type II GRBs are simply associated 
with "long/soft/long lag" GRBs while Type I GRBs are 
simply associated with "short /hard/short lag" GRBs. In 
particular, we will address the nature of the Other SGRB 
Sample. It is likely that some (maybe many) GRBs in 
the Other SGRB Sample are associated with Type I. The 
question is whether they are ALL associated with Type 
I. 

6.1. Are all short/hard GRBs associated with Type I? 

The most straightforward possibility is to accept that 
all short/hard GRBs are associated with Type I GRBs. 
Inspecting the Other SGRB Sample, one may raise 
the following arguments in support of this suggestion: 
(1) They indeed occupy the short/hard domain of the 
BATSE T90-HR diagram, which is in distinct contrast to 
Type II GRBs that predominantly occupy the long/soft 
domain; (2) Most of them deviate from the E p — E^^ so 
relation and the 1? iso — lag relation for Type II GRBs; 
(3) The redshift and luminosity distributions of the ob- 
served sample are different from those of Type II, al- 
though with much overlap; (4) The afterglow luminosi- 
ties are systematically lower than those of the Type II 
majority, although with some overlap. However, as dis- 
cussed below, there are reasons to be suspicious of this 
straightforward conclusion. 

6.2. Are some short/hard GRBs associated with Type 

II? 

The fact that some Type II Gold Sample GRBs are in- 
trinsically short naturally raises the possibility that the 
observed short/hard GRBs are contaminated by Type 
II GRBs. A small contamination is expected given the 
overlapping log-normal distributions of Tqq for the two 
populations. A more intriguing possibility is that the 
contamination is not the simple extension of the Tgo dis- 
tributions, but accounts for a good fraction of the ob- 
served short /hard GRBs. 

Conservatively speaking none of the 4 arguments dis- 
cussed in §5.1 is conclusive. Comparing with the Type 
I Gold Sample, one cannot straightforwardly demon- 
strate that the Other SGRB Sample and the Type I 
Gold Sample come from the same parent sample. In 
particular, the Type I Gold Sample GRBs are relatively 
"long" and "soft" within the short /hard population, with 



4 out of 5 having extended emission. The more "classi- 
cal" short/hard ones, such as GRBs 051221 A, 060313, 
and 061201, do not satisfy the Gold Sample criteria. 
After the 2005 revolution of discovering GRB 050509B 
and GRB 050724, which are associated with elliptical or 
early type host galaxies, it is now clear that such asso- 
ciations are not very common. Most short/hard GRBs 
are found to be associated with star forming galaxies, 
some of which have pr operties close to those of Type 
II GRBs. According to iBergerl (|2009f) . the majority of 
short GRBs appear to occur in star forming galaxies. 
Although some compact merger events can have short 
merger ti me scales, and therefore appear in star forming 
galaxies dBelczvnski et al. 2006; Zheng fc Ramirez-Ruizl 
l2007tlO'Shaughnessv et al.ll2008D . the dominance of star- 
forming host galaxies of short GRBs may raise a con- 
cern regarding whether some short GRBs may be asso- 
ciated with Type II. Although there is an observational 
selection effect that favors redshift identification for star- 
forming galaxies, the GRBs with bright host galaxies 
but no redshift identifications (e.g. GRB 051210, La 
Parola et al. 2006; Berger et al. 2007) are not the 
predominant population of short/hard GRBs. The frac- 
tion of Type II contamination can be small (e.g. < 1/3 
according to Berger 2009), but may not be negligible. 
We speculate that some high-L short GRBs in star- 
forming galaxies may instead be Type II GRBs. Fo r ex- 
ample, GRB 060121 (|de Ugarte Postigo et al.ll2006h has 
E Jtiso = 2.2 x 10 53 erg and L 7iPiiso = 3.4 x 10 54 erg s" 1 
for z = 4.6 and E 1: i SO = 4.2 x 10 52 erg and L 7iPi i SO = 
2.4 x 10 53 erg s^ 1 for z = 1.7. These energy values are 
typical for Type II GRBs, and the luminosity values even 
belong to the bright end of the Type II distribution. Its 
optical afterglow luminosity is also typical for Type II 
and much brighter than those of the Type I Gold sample. 
Given the possible redshifts, this burst lies right on the 
Amati-relation and the luminosity-spectral lag relations 
of most Type II GRBs (see Figs. [3] and @]). GRB 060313 
(|Roming et al.ll2006[) . whose z < 1.1, can have E^^ so ~ 
3.4 x 10 52 erg s _1 for z — 1. Some afterglow light curve 
features (e.g. the very shallow decay of the UVOT, 
Roming et al. 2005, light curve with flickering features) 
are hard to accommodate within the merger scenarios. 
GRB 061201 also s hows a very flat early light curve 
(jStratta et alJ l2007h in both X-ray and optical bands, 
with the first clear X-ra y "plateau" appear ing in a short 
GRB. Some models (e.g. lKumar et al.ll2008f) attribute X- 
ray plateaus to the signature of massive star accretion. 
Within such a scenario, GRB 061201 is then a Type 
II candidate. GRB 060121 al so shows strong tempo- 
ral va riability in the afterglow (jde Ugarte Postigo et al.l 
I2006h . similar to GRB 060313. Even for the not very 
energetic short/hard GRB 051221 A, its identity as a 
Type I GRB is not unquestionable. Its host galaxy 
has log(SSFR)= 0.804, which is greater than those of 
many GRBs in the Type II Gold Sample (Savaglio et al. 
2009). A jet break is detected through ToO observa- 
tions with the Chandra X-ray Telescope (|Burrows et alJ 
l2006h . which gives a jet corrected gamma-ray energy of 
E 1 ~ 5 x 10 49 ergs. This is smaller than but not far off 
from the distribution of Type II GRBs dFrail et al 
200lMBloom et al.ll2003l : iLiang et afl2008tlRacusin et al. 
2009h . and it is also much higher than some other Type 
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I GRBs (e.g. GRB 050509B 
10 48 ergs, Gehrels et al. 2005). 



with E lt i SO r*> 1.1 x 
Although a low den- 



sity n ~ 10 



10 cm was inferred from a fterglow 



modeling (Sod erberg et al.ll2006l : [Burrows et al.ll2006f ). it 
still belon gs to the reasonable n-range of other Type 
II GR Bs |Panaitescu fc Kumarl I2001L 120021 : lYost et all 
120031) . The recent short/hard GRB 090510 detected 



by both Swift and Ferm i (GBM/LAT) fHoyersten et al 
2009t lOhno et all 120091 : iGuiriec et all l2009t iRau et al 



20091 ) is also located in a star forming host galaxy, 



and has an inferred total (collimation-corrected) jet 
gamma-ray and kinetic energies (assuming n ~ 1 cm~ 3 ) 
E 1 /Ek ~ 10 50 erg. It is again not f ar-off from 
the j£y / E K distributions of Type II GRBs dFrail et al 
200lMBloom et al.ll2003t iLiang et al]|2008tlRacusin et al 
and a Type II origin is possible. Finally, 
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Nyse wander et all (|2009l ) pointed out that the optical-to- 



X-ray flux ratios of short GRBs are quite similar to those 
of long GRBs, suggesting a similar circumburst medium 
density for the two populations. This is consistent with 
most short GRBs being associated with Type II. 

Another strong argument for the Type-II association 
of some (even many) short GRBs is related to luminos- 
ity function analyses. Following the similar methodol- 
og y of modeling L — z distribution of Typ e II GRBs 
in IVirgili et all (|2009af ). IVirgili et alJ (|2009bh have stud- 
ied the required luminosity function of Type I GRBs in 
order to reproduce the observed L — z distribution for 
both Type I and Other SGRB samples. The results sug- 
gest that the underlying luminosity function (defined as 
N(L)dL oc L~ q ) must be very shallow (e.g. q ~ 0.5) in 
order to reproduce the L — z distribution data. This shal- 
low luminosity function is different from Type II GRBs 
and other astrophysical objects. A more severe problem 
is that it cannot reproduce the observed log N — log P 
distribution of BATSE short/hard GRBs. This appar- 
ent conflict disfavors the hypothesis that all short/hard 
GRBs are associated with Type I. 

Theoretically, the duration of a GRB is defined by 
Eq.JT]). We now discuss the three relevant time scales 
in turn and address how a short GRB can in principle be 
associated with Type II. 

Firstly, recent studies of the collapsar model suggest 
that the engine time scale t eng i ne may not be al ways long. 
Acco r ding to the standard collap s ar model (Woosley 
119931: iMacFadven fc Wooslevi [l999l : iProga et alJ l2003f ), 
tengine can last as long as the fallback material from the 
collapsar envelope is available to fuel the accretion disk 
or torus. However, one should bear in mind that the 
rotating torus may form only when the specific angu- 
lar momentum of the accreting gas is higher than the 
so-called critical specific angular momentum value, i.e. 
i cr it = 2i? g c, where i? g is the gravitational radius. Note 
that Z cr it is proportional to the mass of the BH. During 
the collapsar evolution the mass accretion rate is very 
high, therefore the BH mass and consequently the crit- 
ical angular momentum increase very fast. As a result, 
the specific angular momentum of the rotating material, 
which was initially sufficient for the torus formation (i.e., 
when the BH was just formed) , may become insufficient 
at a later stage of the collapsar evolut ion when the BH 
mass increases. Uaniuk fc Progal (|2008l ) showed that the 
simple, often cited, estimates of the total mass avail- 



able for torus formation and consequently the duration o f 
a GRB (|MacFadven fc Wooslevi H999I : IProga et al.ll2003h 
are only upper limits. They revised these estimates by 
taking into account the long term effect so that as the 
BH accretes the minimum specific angular momentum 
needed for torus formation increases. These new esti- 
mates predict a significant (an order of magnitude) re- 
duction of the total energy and overall duration of the 
central engine tengine because only a fraction of the ro- 
tating stellar envelope can form a torus. 

If a Type II GRB is po wered by the black hole spin 
(|Blandford fc Znaieldll977f ) rather than accretion, tengine 
of a Type II GRB can be also short, since accretion of ma- 
terials with a very low specific angular momentum would 
slow down the BH and consequently suppress the jet pro- 
duction. The interplay among the BH mass, BH spin 
parameter, and the critical specific angular momentum 
of accreting gas needed fo r the t orus to form have been 
discussed bv I Janiuk et al.l ((2008) . They studied several 
different cases and reached the conclusion that depend- 
ing on the parameter settings, tengine can be as short as 
a second. 

Secondly, the time scale during which a relativistic jet 
is launched it jet) may be in principle shorter than the 
central engine activity time scale (t engine)- There is no 
working baryon-loading model for GRBs, and it is not 
clear how a clean, high entropy outflow is launched. For 
example, if the engine power has several episodes and 
the power in the earlier episodes is not high enough, 
the earlier jet may be choked or be launched but with 
a heavy baryon loading. The ejecta may therefore be- 
come a "dirty" fireball. The GRB episode is then only 
related to the late "clean" fireball phase when baryon 
loading is reduced. 

Thirdly, energy dissipation is needed to convert 
other fo rms (kinetic and magn e tic) of energy to ra- 
diation dRees fc Meszarosl 119921: jMfezaros et alJ 119931 : 
iRees fc Meszaroslll994lThompsonill994D . For a baryonic 
fireball, a steady outflow may not generate significant 
internal, non-thermal emission without internal shocks. 
The energy dissipation time scale (idis) can be smaller 
than tjet- This gives an additional room to reduce the 
duration of Type II GRBs. 

Finally, some other possibilities of producing short 
GRBs from collapsars ha ye been proposed i n the pre- 
Swift era. For example, IZhang et al.l (|2003l) proposed 
that a short/hard pulse of gamma-ray emission may be 
associate d with eruption of the fi reball from the stellar 
envelope. Yamazaki et al. (2004a]) envisioned a geomet- 
ric model to unify long and short GRBs based on a line- 
of-sight effect. The original pictures proposed in these 
papers are no longer supported by the current data, but 
some ideas may be borrowed to associate short GRBs 
with Type II model category. 

6.3. Is there a "Type III" model category? 

The current data do not demand the existence of a 
third type of GRB models to be associated with cos- 
mological GRBs 21 , i.e. those neither associated with 
massive star deaths nor compact star mergers. How- 

21 Again "cosmological GRBs" do not include SGR giant flares 
that are believed to account for a small fraction of short/hard 
GRBs. 
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ever, the possibility is not ruled o ut by the da ta, ei- 
ther. Some "hostless" short GRBs (iBergerlfe OOQ). some 
long GRBs without X-ray afterglows (|Vetere et al.ll2008f ). 
and the SN-less long-duration low-energy GRB 060505 
(iFvnbo et al.l 120061: K)fek et alJ l2]M iThone et all 12001 
iMcBreen et al. 2008; Kami et al. 2003) are oddballs that 
may hold the clues to identify possible new model cate- 
gories of GRBs. 

7. NATURE OF GRB 080913 AND GRB 090423 

We now discuss the possible origin of GRB 080913 at 
z = 6.7 and GRB 090423 at z = 8.3. We mainly focus on 
GRB 080913. The case of GRB 090423 is amazingly sim- 
ilar to GRB 080913, and the conclusion for GRB 080913 
can be directly applied to GRB 090423 as well. 

7.1. Prompt properties and empirical correlations 

As discussed in §2 (FigOJ, GRB 080913 and GRB 
090423 appear as long GRBs in the observer's frame, 
but are intrinsically short/hard GRBs in the rest frame. 
If one applies the criterion for the observed Tgo, both 
bursts are "long" and therefore may be associated with 
Type II according to the ansatz Eq. (fT7|) . However, the 
association of a particular GRB to a particular physical 
model type should not have a z-dependence. The identi- 
cal burst, if it have occured at z < 1, would be recognized 
as a short/hard GRB, and hence, a Type I candidate ac- 
cording to the ansatz Eq. (fT7|) . So it is not straightfor- 
ward to determine the physical model category a GRB 
should be associated with based on the observed Tgo-HR 
data. 

We inspect the compliance of GRB 080913 and GRB 
090423 with the empirical correlations. First, they are 
consistent with the E p — E 7yiso Amati- relation (FigEJ) , 
although GRB 080913 is near the 3a upper boundary. 
This has been regarded as one arg ument in support of th e 
Type II origin of GRB 080913 bv lGreiner etatl (|2009aD . 
On the other hand, GRB 080913 is also consistent with 
the new track defined by Type I and short/hard GRBs 
within 3cr. This suggests that the possibility that it is 
a Type I GRB (or at least similar to other short/hard 
GRBs) is not ruled out based on this criterion. The com- 
pliance of GRB 090423 with the Amati-relation is more 
robust. Second, inspecting the iso — lag correlation, it 
seems that both GRB 080913 and' GRB 090423 are con- 
sistent with being Typ e II GRBs - this was another argu- 
ment bv lGreiner et al.l (|2009af ). However, both GRBs are 
also consistent with the "zero-lag" trend of short GRBs. 

In conclusion, based on statistical properties, both 
GRBs can be taken as good Type II candidates. How- 
ever, the criterion based on empirical correlations is not 
robust enough to claim the case, and supports from other 
criteria (see below) are needed to draw firmer conclu- 
sions. 

7.2. Afterglow properties 

The rest-frame broadband (X-ray and optical) after- 
glow luminosities of GRB 080913 are moderate, brack- 
eted between those of Type II and Type I Gold Sampels 
(Fig[5] and FigJT]). Although the early luminosities are 
relatively low, a distinct energy injection episode raises 
the afterglow luminosit y level of this burst to those of 
Type II at later epochs. iGreiner et al.l (|2009al ) have mod- 
eled the afterglow and suggested that the light curves are 



consistent with the deceleration of a relativistic jet by a 
dense circumburst medium with constant density. The 
data are consistent with the existence of an achromatic 
plateau in X-ray and several optical/IR bands, which can 
be interpreted within the framework of a continuously 
fed forward shock. Alternatively, the X-ray rebrighten- 
ing around 10 5 s may be due to an X-ray flare, whose 
softer emission may also account for the rebrightcning in 
the IR/optical bands. 

Improving upon IGreiner et al.l (|2009af) , we have per- 
formed a more detailed numerical modeling of the broad- 
band afterglow data. The optical emission can be well 
modeled by standard synchrotron radiation from the for- 
ward shock, while the X-ray emission is likely dominated 
by the synchrotron self-Compton emission (SSC). The 
following parameters can fit the optical data well (al- 
though we do not apply a parameter search to judge 
whether this is the best fit): the initial isotropic kinetic 
energy of the fireball Eic,iso ~ 3.7 x 10 52 erg, the am- 
bient density n ~ 3000 cm -3 , the electron equiparti- 
tion parameter e e ~ 0.04, the magnetic field equipar- 
tition parameter eb ~ 10~ 5 , and the electron spectral 
index p ~ 2.2. If the late rebrightening is interpreted 
as an energy injection from the central engine with a 
time-dependent luminosity L = Lo(t/to) q , the data are 
consistent with to ~ 6.5 x 10 3 s, L ~ 2.9 x 10 50 erg s~ x , 
q ~ 1 for t < to, and q < — 1 for t > t<j. In any case, 
a high-density constant medium is nee ded. This is con- 
sistent with the expectation at high-z (|Gou et al.ll2004f ) 
as well as the fitted density of GR B 050904 at z = 6.3 
(|Frail et all 120061 : IGou et all 12001 . The data also de- 
mand a jet opening angle 9j > 0.22 rad, which corre- 
sponds to a geometrically-corrected total gamma-ray en- 
ergy E 1 > 1.7 x 10 51 ergs, and a geometrically-corrected 
total kinetic energy Ek > 9.0 x 10 50 ergs. This value 
is consistent with those of Type II GRBs dFrail et al. 
200lMBloom et alj|2003t iLiang et alj2008tlRacusin et al. 
2009). This is probably the strongest argument in favor 
of associating the burst with Type II. 

GRB 090423 has even brighter X-ray and optical af- 
terglow luminosities than GRB 080913. Although we 
did not perform detailed afterglow modeling, the after- 
glow parameters favor those of Type II GRBs, similar to 
GRB 080913. 

7.3. Short Type II or high-z Type I? 

Since the discovery of GRB 080913, there has 
been a debate about its progenitor. As discussed 
above, data analyses and theoretical modeling sug- 
gest that GRB 0809 13 is very likely a Type II GRB 
(IGreiner et al.l l2009al) , alth o ugh a Type I a s sociat ion 
(jPerez-Ramirez et al.l 120081: iBelczvnski et al.l I2008D is 
not ruled out. The evidence in support of the Type 
II origin of GRB 080913 includes: large values of 
the geometrically-corrected gamma-ray (E 7 ) and kinetic 
(Ek) energies, moderately bright intrinsic afterglow lu- 
minosities, a required high density of the circumburst 
medium, and the marginal compliance of the E p — E Jy i SO 
relation of Type II GRBs. 

On the other hand, if GRB 080913 were a high-z 
Type I GRB, as suggested by its intrinsically short du- 
ration, it would have to be an energetic merger event, 
likely due to a BH-NS merger with a rapidly rotat- 
ing massive BH. The energy tapping mechanism would 
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have to be the [Blandford fc Znaiekl (|1977f ) mechanism 
(jPerez- Ramirez et al.l 120081 ). For this possibility, one 
requires that during the short age of the Universe at 
z = 6.7, i.e. r ~ 8.3 x 10 8 yr for the concor- 
dance universe, a BH-N S system is formed and merged. 
iBelczvnski et al.l (|2008l ) have modeled this possibility in 
detail, and claimed that the event rates for massive star 
core collapses that give rise to Type II GRBs and for com- 
pact star mergers (both NS-NS and BH-NS) are compa- 
rable at z = 6.7. They concluded that both scenarios are 
possible. However, there are se veral factors that would 
change this conclusion. First, IBelczvnski et al.l (2008) 
assumed that all the mergers give rise to GRBs. In real- 
ity it may be that only a fraction of mergers give rise to 
GRBs. This fraction factor may be calibrated through 
confronting the observed number ratio of Type I and 
Type II GRBs with the model predictions. In the cur- 
rent population synthesis models, this factor is not taken 
into account (K. Belczynski, 2008, personal communica- 
tion). Secondly, GRB 080913 would be a high-luminosity 
Type I GRB if it is associated with that category. Con- 
sidering the p ower law luminosity function inferred for 
Type I GRBs (|Virgili et al.ll2009bf ). detecting one high- 
L event would demand many more low-L events, which 
would require a significant increase of the required event 
rate of compact star mergers that is inconsistent with 
the results of population synthesis. Finally, as we ar- 
gued above, the large value of the geometrically-corrected 
gamma-ray and afterglow energies do not favor a NS-NS 
merger model. Only BH-NS mergers with a highly spin- 
ning BH should be counted. This would greatly reduce 
the theoretically predicted event rate that satisfies the 
constraint. 

The detection of GRB 090423, another intrinsically 
short/hard, high-z, high-L GRB, strongly supports a 
Type II association of both GRB 080913 and GRB 
090423. As argued above, the probability of detecting 
a high-L Type I event at high-z is much smaller than 
that of detecting a moderate- L Type II event. With one 
event, one may still argue for a chance coincidence. With 
the detection of GRB 090423, the chance probability of 
detecting two high-L merger events at high-z is greatly 
reduced, and one can more firmly associate both GRBs 
with Type II. 

With two intrinsically short high-z Type II GRBs de- 
tected, one must ask why these events tend to exist 
at high-z. One possibility would be that it is simply 
a threshold selection effect. Both events have moder- 
ate gamma-ray luminosities, and were detected not far 
above the threshold. It is possible that there exists other 
softer pulses that are below the sensitivity threshold of 
Swift/BAT (J. S. Bloom, 2009, private communication). 
On the other hand, such softer emission would be easily 
detected by Swift/XRT if it was indeed there. Both GRB 
080913 and GRB 090423 have X-ray flares. However, ex- 
trapolating them into the gamma-ray band using a sim- 
ple spectral model suggests that they would appear as 
low- level extended emission of short GRBs (Fig.l). One 
may still argue for missing soft emission before the XRT 
slew. This is not ruled out, but the XRT slew time cor- 
responds to a rest frame time 99.5/7.7=12.9 s for GRB 
080913 and 72.5/9.3=7.8 s for GRB 090423. The intrin- 
sic Tgo's have to be in any case smaller than these values 
for these bursts. 



A more intriguing possibility would be that this is due 
to a physical origin and reflects the intrinsic property of 
the high-z massive stars. These high-z stars may not be 
spinning as rapidly as their low-z sisters, so that only a 
smaller mass is left after the prompt collapse. More high- 
z GRB data are needed to test whether such a scenario 
is demanded by the data. 

8. HOW TO ASSOCIATE A BURST WITH A PHYSICAL 
MODEL CATEGORY? 

The extensive discussion presented above suggests that 
it is not always easy to associate a particular GRB to 
a particular physical model category based on observa- 
tional criteria. The multiple observational criteria dis- 
cussed in this paper are su mmarized in T able 2. This is 
an extension of Figure 2 of I Zhang] f2006). New criteria 
are added based on the discussion in this paper. A new 
column lays out the issues of each criterion. The criteria 
are sorted by relevant observations. The first six rows 
(duration, spectrum, spectral lag, Ey t i so , E p — E~ ( ^ so re- 
lation, and L^ iso -lag relation, are based on the gamma- 
ray properties only. The next five rows (supernova as- 
sociation, circumburst medium type, Ex,iso, jet opening 
angle, and the geometrically corrected energies L 7 and 
Ek ), are based on follow-up broadband observations and 
afterglow modeling. The next three rows (host galaxy 
type, specific star forming rate of the host galaxy, and 
offset of the GRB from the host galaxy) are based on 
observations of the host galaxies. The next two rows 
(redshift distribution and luminosity function) are statis- 
tical properties. The final row is the gravitational wave 
criterion. In general, most of these criteria are not "con- 
clusive", i.e., one cannot draw a firm conclusion based 
on a single criterion. Nonetheless, there are several cri- 
teria which, if satisfied, would unambiguously associate 
a GRB to a certain physical model category. These are 
marked in bold in Table 2. In particular, if a GRB is 
found in an elliptical or an early type galaxy, or if the 
SSFR of its host galaxy is very low, one would be able 
to associate it with Type I. On the other hand, a SN 
association or the identification of a wind-type medium 
in a GRB would establish its association with Type II. 

Unfortunately, the above four criteria are usually not 
satisfied for most GRBs. One is then obliged to use mul- 
tiple criteria since there are overlapping predicted prop- 
erties between the two physical model types for each in- 
dividual criterion. In Fig[8] we cautiously propose an 
operational procedure to discern the physical origin of a 
GRB based on the available data. 

Several features are worth commenting on in Fig[8] (1) 
The criteria to define the physical model category a burst 
is associated with are less stringent compared with those 
used to define the Gold Samples in §5.1. This is because 
the purpose of the Gold Samples was to allow us to per- 
form statistical analyses. After reviewing the statistical 
properties in §5, we have gained confidence on additional 
criteria so that more bursts can be analyzed. (2) There 
are five outcomes in the flowchart. Besides the solid Type 
I/II identifications, we also define Type I/II "candidates" 
and the "unknown" category. The Type I/II candidates 
refer to those with evidence of associating a burst to a 
particular physical model category, but the evidence is 
not strong enough to make a firm claim. The unknown 
category includes the oddball GRBs that do not obvi- 
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ously fit into any criteria discussed in this paper, or the 
observational data are not adequate for us to make the 
judgement. They may be associated with Type I, Type 
II or a completely new type of models. (3) Some qualita- 
tive rather than quantitative criteria have been used (e.g. 
high/low SSFR, large offset, large/small E 1 , E K )- The 
reason is that it is very difficult to adopt quantitative 
criteria at the current stage, since the distributions of 
these quantities predicted by both physical model types 
and displayed in the statistical analyses of the Type I/II 
Gold Samples are continuous, without sharp transitions. 
The "high/low" and "large/small" definitions are based 
on the statistical properties, and therefore in the relative 
sense. If confusion occurs (e.g. the quantity is near the 
boundary and not easy to judge whether it is high/low, 
large/small, one can follow the "?" sign to go down 
the flowchart. The flowchart is reasonably operational, 
i.e. essentially every GRB with reasonable afterglow 
follow up observations can find a destiny in the chart. 
For example, the SN-less long-duration GRB 060614 
(iGehrels et al.ll200& iGal-Yam et all I20M lFvnbo"etal] 
120061: iDella Valle et aljl2006al ) is associated with Type I 
(based on low SSFR), and the other SN-less GRB 060505 
(iFvnbo et alJl20M lOfek et al.ll2007t iThone et all 12003 
iMcBreen et all |2008|T can be associated with a Type I 
candidate based on its small energetics, or an "unknown" 
burst if one argues that the iso — lag relation is satis- 
fied for this burst (jMcBreen et al.l I2008D . GRB 080913 
and GRB 090423 (the main topic of this paper) find their 
homes as Type II candidates based on the E p — Ej^ so 
correlation. GRB 060I2I (a high-z short GRB) satisfy- 
ing the E p — E~ t i so is also found to be associated with 
the "Type II candidate" outcome in the flowchart. (4) 
It is possible that the procedure and the criteria may be 
further revised as more data are accumulated. The cur- 
rent procedure only reflects the best knowledge for the 
time being. 

In the flowchart, there are five thick arrows that bridge 
the short-duration and long-duration GRBs. This sug- 
gests that the duration information sometimes is mislead- 
ing. Some long duration GRBs can be associated with 
Type I (e.g. GRB 060614 and probably GRB 080503, 
Peyley et al. 2008), and some short duration GRBs can 
be associated with Type II (e.g. GRB 060121, GRB 
080913 and GRB 090423). We also present two dashed 
arrows in the flowchart. These two tracks (a short GRB 
associated with a SN and a long GRB with an ellipti- 
cal/early type host galaxy) are in principle possible, but 
such bursts have never been observed so far 22 . The or- 
der of the criteria in Fig. 8 is based on the "defmiteness" 
of the criteria, with the higher-level ones carrying more 
weight than the lower-level ones. Notice that "hardness" 
is generally not regarded as a definitive criterion in the 
flowchart (except for the relative hardness of the short 
spike and the extended emission). 

9. SUMMARY 

Prompted by the interesting question whether the 
z = 6.7 GRB 080913 and z = 8.3 GRB 090423 are 

22 The GRB field is full of surprises. If some short/hard GRBs 
are indeed associated with Type II as argued in this paper, one 
may someday discover a SN associated with a short/hard GRB. 
We encourage continuous SN searches for all nearby GRBs, both 
long and short. 



intrinsically short GRBs associated with th e Type II 
physical model category (jGreiner et alJl200 9a) or high-z 
GRBs associated with the Type I physical model cate- 
gory (|Perez-Ramirez et alj|200a) . we performed a more 
thorough investigation on the two physically distinct cat- 
egories of GRB models and their predicted observational 
characteristics. We further developed the "Type I/II" 
concept proposed in lZhang et alj (|2007bl ) in the following 
directions. (1) We have reviewed and expanded the pos- 
sible multiple observational criteria, and discussed their 
physical origins from the theoretical point of view. By 
doing so, we are able to differentiate those criteria that 
are more closely related to the progenitor types and those 
that are more directly related to radiation physics. In 
particular, we argue that SN association, host galaxy 
properties (type and SSFR), and the offset of the GRB 
location in the host galaxy are more directly related to 
the progenitor types. The gamma-ray properties, such 
as duration, hardness, spectral lag, empirical correla- 
tions, are more related to jet dissipation and radiation 
processes in the emission region, and can only be re- 
lated to progenitors indirectly. Afterglow and statisti- 
cal properties can be used to diagnose GRB progenitor, 
but theoretical modeling is needed. Gravitational wave 
signals may be the best criterion to directly probe the 
progenitor system, but they are too faint for the cur- 
rent detectors to detect. (2) We use several key observa- 
tional criteria that are directly related to GRB progeni- 
tors to define the Gold Samples for Type I and Type II, 
respectively. These criteria do not involve GRB gamma- 
ray emission properties such as duration, hardness, spec- 
tral lag, etc. We then use these samples to investigate 
their statistical properties, especially their distribution 
in the duration-hardness space. We found that the Type 
II Gold Sample represent the BATSE long/soft popula- 
tion well. The Type I Gold Sample, on the other hand, 
is not very representative of the short/hard population. 
The Type I Gold Sample GRBs are typically "long" and 
not particularly "hard". (3) Although some short/hard 
GRBs detected in the Swift era may share a similar ori- 
gin as the Type I Gold Sample, we suggest that some 
(maybe most) high-L short GRBs may be instead asso- 
ciated with Type II, namely, of a massive star origin. 
(3) We summarized the multiple observational criteria 
needed to discern the physical origin of a GRB in Table 
2, with various issues laid out. We emphasize that it is 
not always straightforward to judge the physical model 
category a particular GRB is associated with, and we 
cautiously proposed an operational procedure to discern 
the physical origin of GRBs (FiglHJ) - (4) According to 
this procedure, GRB 080913 and GRB 090423 are Type 
II candidates. Although a specific Type I scenario invok- 
ing the Blandford-Znajek mechanism of a BH-NS merger 
system is not completely ruled out, the fact that two such 
GRBs are detected at high-z indeed suggest that a Type 
I association of these bursts is essentially impossible. 

The proposed procedure to associate a particular GRB 
to a particular physical model category is subject to fur- 
ther test with new observational data. More detailed 
analyses may allow more quantative criteria to discern 
the physical origin of GRBs. Based on past experience, 
the chances are high that new observations will bring 
surprises that continuously call for modifications of the 
criteria, which would further our understanding of the 
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physical origins of cosmological GRBs. 
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TABLE 2 

Observational criteria for physically classifying GRBs. 
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Fig. 8. — A recommended procedure to judge the association of a particular GRB to a particular physical model category. Multiple 
observational criteria have been applied. Question marks stand for no information being available to judge the validity of the criterion. 
The two dotted arrows stand for the possibilities that are in principle possible but have never been observed. Five thick arrows bridge the 
long-duration and short-duration GRBs, suggesting that the there can be long duration Type I and short duration Type II GRBs. 



APPENDIX 
DETAILS OF THE TYPE I GOLD SAMPLE 

• GRB 050509B: No optical afterglow is detected. The host galaxy is very likely a bright cD elliptical galaxy in a 
nearby galaxy cluster at z = 0.2248 ±0.0 002. The Swift XRT error circle is offset from the elliptical ga laxy. Deep 
upper limit on SN association achieved (|Gehrels et al. 1 120051 iBloom et al.l 120061 ; iHiorth et al.l l2005bf) . Satisfies 
both criteria # 1 and 2 of the Gold Sample. 

• GRB 050709: The host galaxy is identified with the optical afterglow, and is a star forming galaxy at z = 
0.1606 ± 0.0001. The loc ation of the GRB is at the outskirt of the host. Deep upper limit on SN association 
achieved (|Fox et al.ll2005| ). Satisfies criterion # 2 for the Gold Sample. 



• GRB 050724: The host is an early type galaxy at z = 0.2576 ± 0.0004. The afterglow is at the outskirt of the 
host (Barthcl mv et al.ll2005aHBerger et~a l. 2005). Deep upper limit on SN association achieved ( |Malesani et al.l 
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l2007t) . Satisfies both criteria #1 and 2 for Gold Sample. 

• GRB 060614: The host galaxy has a low SSFR (|Gal-Yam et all 120061: ISavaglio et all I2009D. Very stringent 
upper limits for any ass ociated SN fainter than any known SN achieved (jGal-Yam et al.ll2006l : iFvnbo ^eTaTII^OOa 
Deli a Valle et al.ll2006al ). Satisfies criterion # 2 for the Gold Sample. 

• GRB 061006: A faint host galaxy is det ected at z = 0.4377 ± 0.0002 (jBerger et al.ll2007t) . The SSFR is very low 
(jBerger et alj|200"7t Savaglio ct al. 2009). Satisfies criterion # 1 for the Gold Sample. 

DETAILS OF OTHER SHORT/HARD GRB SAMPLE 

• GRB 000607: This was a IPN-localized GRB. A putative host galaxy at z — 0.14 was proposed (|Gal-Yam et al.l 
l200l . 

• GRB 050 813: No optical after glow detected. The X-ray error circle is associated with a galaxy cluster at 
z ~ 0.72 (jProchaska et al.|[2006f h Moderately deep SN limit (relevant to this low-z interpretation) was reported 
by (jFerrero et alll2007ft . 

• GRB 051210: A Swift GRB with X-ray afterglow (|La Parola et al.|[2006T ). A galaxy appeared outside the error 
box, but is likely the host. No lines are observed. It is argued that z > 1.4 (jBerger 

• GRB 051221A: Host gala xy is a star forming galaxy at z = 0.5464 with a very high SSFR (jSoderberg et alj 
[2006l :[S avaglio et al. 1 12009). A bright S N such as 1998bw is ru led out, but the limit is still consistent with the 
existence of faint SN 2002ap-like event (jSoderberg et al.l l2006). So it does not satisfy the citeria of the Type I 
Gold Sample. 

• GRB 060121: HETE-2 GRB with a faint optical afterglow, leading to the di scovery of an extremely faint 
host galaxy The redshift of the afterglow can be estimated as either 4.6 or 1.7 (|de Ugarte Postigo et alj|2006t 
iBerger et aT1l2007| ). 

• GRB 060313: Bright short/hard GRB with bright af terglow (|R oming et al J 12006^ 1 . A faint host galaxy is iden- 
tified whose redshift is unknwon (jBerger et alJ 12007). Spectral analysis of the UVOT data suggests z < 1.1 
(jRoming et al.ll2006h . 

• GRB 0605 02B: No optical af terglow is detected. The XRT position is close to a nearby early type galaxy at 
z = 0.287 (HI 00m et aT]|2007f ). The chance probability for the associa tion is 0.03. There is another faint object 
in the field of view, which could be the host galaxy at a high redshift (jBerger et al.ll2007t) . 

• GRB 060801: No optical afterglow is detected. Two possible sources may be considered as the host. One has a 
redshift z = 1.131 (which is slightly outside the XRT error box using the UVOT-a ligned XRT position). Another 
source is within the error box, but is likely even farther away (jBerger etaD[2QQ3). 

• GRB 061201: Optical afterglow was detected by UVOT. No host galaxy was identified. Candidates include 
galaxy cluster Abel l 995 (z = 0.0865), a star-forming galaxy at z — 0.111, or a missing host at an even higher z 
(jStratta et aT]|2007ft . 

• GRB 061210: A Swift GRB with delayed X-ray afterglow detection. No optic al afterglow detected. The host 
galaxy is likely a star forming galaxy at z = 0.4095 ± 0.001 (jBerger et al.ll2007T ) ■ 

• GRB 061217: A faint Swift burst without optical aft erglow detection. Within the XRT error circle, there is a 
star forming galaxy at z = 0.8270 (jBerger et al.ll2007| ). 

• GRB 070429B: A Sw ift GRB with delayed X-ray afterglow d etection. The host galaxy is likely a faint galaxy at 
z = 0.9023 ± 0.0002 (jPerlev et alJl2ot)7al : iCenko et al.ll2008al K 

• GRB 070707: Detected by INTEGRAL and have X-r ay and optical afterglow detected. A very faint host galaxy 
candidate was reported with no redshift information (jPiranomonte et aMl2008f ) . 

• GRB 070714B: A Swift GRB with opti cal afterglow. A secure host galaxy at z = 0.9225 ± 0.0001 is identified 
(jGraham et alll2009t [Cenko et al.ll2008al ). 

• GRB 070724A: A S wift GRB with X-ray a fterglow. A potential host galaxy is detected, which is a star forming 
galaxy at z = 0.457 (jCucchiara et al.ll2007T) . 

• GRB 070729: A Swift GRB with faint X-ray afterglow. A putative red host galaxy is identified (jBerger fc Murravl 
I2007D . No redshift is known. 

• GRB 070809: A Swift GRB with X-r ay and optical afterglow. A nearby, edge-on spiral galaxy may be the host, 
with z = 0.2187 (jPerlev et al.ll2007bl ). 
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• GRB 071227: A Swift GRB with X-ray and optical afterglows. A hos t galaxy is identified as an edge on spiral 
galaxy. The redshift is z = 0.3940 (|Bergeril200a iD'Avanzo et al.ll2009ft . 

• GRB 080123: A Swift GRB with X-ray afterglow. No host galaxy detection is reported. 

• GRB 080503: A Swift GRB with short initial spike and very bright extended emission. X-ray and optical 
afterglows are detected. There is no galaxy directly at the GRB position. There are faint gala xies nearby, but 
one cannot make firm statements regarding their association with the GRB (|Perlev et al.ll2009ft . 

After this work is finished, two more interesting short/hard GRBs were detected, whose observational properties 
strengthen the main theme of this paper. We include them here as follows. 

• GRB 090426: This is rest-frame 0.35 s short/hard GRB at z = 2.6 (jLevesque et al.ll2009D . It has a blue, very 
luminous, star-forming putative host galaxy with a small angu lar offset of the afterg low location from the center. 
It is very likely associated to Type II as argued in this paper (|Levesaue et al.l 12009). 

• GRB 090510: This is br i ght short/hard GRB detected by both Swift and Fermi (GBM/L AT) dHoversten et all 
20091: lOhno et al.l 120091: iGuiriec et"aLl I2009D with bright X-ray a nd optical after glows (jGrupe et al.l 12009b: 
Kuin et all [2009t lOlivares et all 12009ft . With a redshift z ~ 0.9 (|Rau et all 12009). the isotropic energy and 
luminosity of this burst all belong to the high end of the distribution for Other Short/Hard GRBs presented in 
Fig. 5. The afterglow is consistent with a uniform density medium. The X-ray lightcurve shows an early break 
at t ~ 1500 s since the trigger with a post-break decay index ~ —2.16. If it is interpreted as a jet break, then 
the total (collimation-corrected) jet gamma-ray and kinetic energies (assuming n ~ 1 cm~ 3 ) are E 1 /Ek ~ 10 50 
erg. This is relat i vely small as compared with (but not far-off from) the E ^/Ek distributions of Type II GRBs 
(jFrail et al.ll2001l : lBloom et alJl200a iLiang et al]l2008t iRacusin et al.ll2009ft . Although one may argue for a Type 
I association based on this, a Type II assoication is not strongly disfavored. 



